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CORRIGENDA  AND  ADDENDA 


vi  and  vil 


add  "in  the  moving  frame"  after  the  word  "frequency"  for  the 
notations  of  ft  and  w^., 

Eq.  (2.1.8)  -  last  term  -  "  5  ^tg  -  t  +  J=Lz_=JL^  " 
should  read  6  ^  tg  -  t  +  I— 1  ^ 


Eqs.  (2.2.8),  (2.2.10),  and  (2.2.12)  -  "(Txx)2  " 

should  read  "  T  ^  " 
xx 

second  last  line  -  "...Davies  (Ref.  36).  This..." 

should  read  "...Davies  (Ref.  36),  this..." 

second  last  paragraph  -  line  3  *  the  word" derivations" 
should  read  "deviations" 

00 

r 

second  equation  of  (7.l)  -  there  should  be  a  "  I  "  sigr 
between  — ^  and  seR(j^,4mi  +  tJcT,0,0,i)  0 


Pig.  2 


last  line  -  "...90°.  The  acoustic..."  should  read 
"...90°,  the  acoustic..." 

add  0-5  after  "AMPLIFICATION" 


Appendix  B  second  paragraph  -  line  1  -  "stationary  in  time,"  sho;ild  read 

"...statistically  stationary", 

Some  additional  remarks  appear  necessary  regarding  the  justifi¬ 
cation  in  extrapolating  our  data  for  the  cases  of  ' shear’  noise  mentioned  or 
page  18.  We  have  pointed  out  there  that  the  fourth  derivatives  at  t  -  0  were 
quite  sensitive  to  the  range  of  data  available.  The  reason  being  that  the 
fitted  curves  behaved  properly  only  within  the  range  of  data  being  fitted.  The; 
would  have  maxima  and  minima  outside  this  range  because  of  the  inherent  nature 
a  polynomial  with  a  mixture  of  positive  and  negative  coefficients.  Ihysicaliy 
speaking  the  curves  in  Fig.  29  should  not  have  these  sinuosities  at  large  r  *s 
indicated  by  the  relatively  smooth  'self'  noise  data  in  Fig.  30.  The  purpose  o 
the  extrapolation,  then,  is  to  ensure  that  the  function  is  of  a  physically  ad¬ 
missible  shape  in  the  extrapolation  range.  The  analytical  function  obtained  by 
least  squares  to  fit  the  extrapolated  data  fits  the  original  data  points  as  wel 
as  before,  and  in  addition  it  now  shows  a  proper  behavior  at  larger  values  of  t 
argument . 


With  the  'shear  ’  noise  curves  so  extrapolated,  the  data  wege 
thought  to  be  on  a  comparable  basis  with  the  'self  '  noise  data;  that  is,  the 
abscissa  scale  is  long  enough  to  bring  the  curves  (which  are  integrals,  of 
correlations)  down  by  about  half.  Comparable  accuracy  of  the  fourth  derivative 
is  then  anticipated. 


It  should  be  pointed  out  here  that  although  the  absolute 
magnitude  of  the  fourth  derivatives  at  t  =  0  and  their  'half  wave-length' 
were  different,  the  relative  magnitudes  between  the  three  cases  of  'shear' 
noise  and  the  general  shape  of  the  fourth  derivative  curves  remained  the 
same  whether  we  have  extrapolated  our  data  or  not. 

Finally  it  is  to  be  understood  that  the  significance  of 
the  inferred  result  at  90  in  Fig.  36  depends  on  the  accuracy  in  obtain¬ 
ing  the  relative  magnitude  between  the  'shear'  and  the  ’self’  noise. 
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SUMMARY 


-  Lighthill’s  equation  for  aerodynamic  noise  has  been  reformu¬ 
lated  in  terms  of  its  spectral  characteristics  using  the  one- dimensional 
fourier  cosine  transform.  The  final  formalism  is  more  revealing  in  that  both 
the  Doppler  effect  of  moving  eddies  and  the  extent  to  which  retarded  time 
differences  can  be  neglected  are  explicit  in  the  final  equation.  Also  the 
convection  effect  can  be  distinguished  as  a  combined  effect  of  Doppler  shift 
and  amplification. 

Experimentally,  two-point  space-time  correlations  of  both 
the  turbulent  velocities  and  the  square  of  these  velocities  have  been  measured 
in  tne  mixing  region  ( rj-^/D  =  4,  t)2/D  =  0.5)  of  a  4  inch  model  jet  (about  142 
fps)  with  constant  temperature  hot-wire  anemometers.  These  measurements  in¬ 
cluded  the  ordinary  u-component  and  also  the  components  at  45°  and  60°  to 
the  jet  axis.  Results  for  the  u-component  agree  well  with  those  of  previous 
investigators . 

Using  the  moving  axis  transformation  of  these  space-time 
correlation  functions  to-gether  with  the  mean  velocity  measurements,  an  esti¬ 
mate  of  the  basic  directivity,  the  intensity,  and  the  spectrum  of  both  the 
’shear'  and  the  'self'  noise  generated  by  a  unit  volume  of  jet  turbulence  was 
obtained.  The  calculations  were  based  on  the  Proudman  form  of  Lighthill's 
equation  and  Ribner's  notion  of  'shear'  and  'self'  noise.  These  results 
show  reasonable  quantitative  comparison  with  actual  acoustic  measurements 
of  other  investigators  and  they  are  also  in  qualitative  agreement  with  some 
of  the  theoretical  predictions  made  by  Ribner. 
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I.  INTRODUCTION 


A  problem  in  connection  with  turbulence  which  is  of  considerable 
interest  and  importance  is  that  of  determining  the  acoustic  radiation  from  the 
turbulent  motion  especially  from  turbulent  jets.  This  'aerodynamic  noise'  pro¬ 
blem  was  first  formulated  successfully  by  Lighthill  (Refs.  1  and  2)  using  an 
acoustic  analogy.  His  work  stimulated  a  great  number  of  workers  in  the  field 
to  extend,  re-examine  and  apply  the  theory.  An  account  of  the  theory  which 
gives  perhaps  the  broadest  coverage  of  the  significant  extensions  applicable 
to  Jet  noise  is  provided  in  a  review  article  by  Ribner  (Ref.  5)*  Comprehen¬ 
sive  surveys  of  experimental  results  are  given  in  Refs.  4,  5>  and  6.  Most  of 
these  previous  experiments  deal  with  the  actual  acoustic  measurements  both  of 
model  jets  and  jet  engines.  There  are  relatively  few  measurements  of  the 
turbulent  flow  in  model  jets  that  are  applicable  to  noise  estimation  (Refs.  7» 

8,  9>  and  10).  Even  these  experiments  are  of  comparatively  simple  types. 
Although  they  have  contributed  a  great  deal  to  the  understanding  of  the  tur¬ 
bulent  structure  of  jet  flows,  they  are  inadequate  for  determining  the  acoustic 
radiation  without  relying  on  major  simplifying  assumptions  to  the  original 
theory. 


1.1  Review  of  Previous  Work 

In  order  to  predict  the  acoustic  radiation  from  known  properties 
of  the  turbulence,  one  has  to  measure  or  deduce  the  two-point  space-tine 
correlation  of  the  double  time  derivative  of  the  Reynolds  stresses  and  inte¬ 
grate  over  the  turbulent  field.  A  rigorious  evaluation  of  this  kind  in  three 
dimensions  would  be  extremely  difficult  if  not  impossible.  The  only  circum¬ 
stance  in  which  the  noise  emission  can  be  explicitly  calculated  is  for  the  case 
of  a  decaying  isotropic  turbulent  field,  considered  by  Proudman  (Ref.  11). 

But  this  particular  case  has  no  direct  application  to  jet  noise  estimation. 

So  far,  there  are  two  comprehensive  developments  of  approximation  to  Lighthill' s 
equation  that  allow  the  noise  emission  to  be  computed  in  terms  of  simple  pro¬ 
perties  of  a  turbulent  jet  such  as  turbulence  intensities,  length  scales, 
moving- frame  time  scales  and  mean  shear.  These  are  the  works  of  Lilley  (Ref. 

12)  and  Ribner  (Ref.  13).  However,  both  treatments  rely  on  the  postulation 
of  homogeneous  and  isotropic  turbulence  and  also  on  the  assumption  that  the 
space-time  correlation  function  of  the  Reynolds  stresses  are  separable  in 
space  and  time. 


Lilley' s  approach  is  based  on  an  approximation  of  Lighthill 's 
equation  and  is  applicable  only  in  regions  of  strong  shear.  For  regions  with¬ 
out  strong  shear  he  used  Proudman' s  result  and  labelled  the  two  types  of  noise 
the  'shear'  and  'self'  noise.  With  Laurence's  results  (Ref.  7)»  he  is  able  to 
show  that  the  noise  sources  are  concentrated  in  the  middle  of  the  mixing  layer 
and  that  slices  of  jet  within  four  diameters  of  the  nozzle  exit  emit  the  same 
noise  power j  beyond  that  the  emission  falls  off  very  abruptly.  This  is  in 
agreement  with  the  theoretical  investigation  by  Ribner  (Ref.  l4)  and  a  deduct¬ 
ion  by  Dyer  (Ref.  15)  based  on  acoustic  measurements. 

Ribner 's  calculation  (Ref.  13)  starts  with  the  Proudman  form 
of  Lighthill' s  equation  which  is  exact.  He  also  treated  the  total  noise  as  a 
sum  of  the' shear'  and  'self'  noise  by  splitting  the  instantaneous  velocity  into 
a  mean  velocity  plus  a  perturbation.  But  the  treatment  is  coherent  in  itself. 
With  Laurence's  results,  he  is  able  to  obtain  a  basic  directivity  of  jet  noise 
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(l  +  cos^S)  that  agrees  better  with  experimental  results  than  the  ' four-leaf - 
clover'  type  (l  +  constant  x  sin226)  as  suggested  by  Lighthill  and  Lilley. 

(See  Figure  1  for  the  definition  of  6.)  Furthermore,  his  treatment  shows  that 
the  'self*  noise  spectrum  peaks  an  octave  higher  than  that  of  the  'shear' 
noise.  With  these  two  spectra,  he  is  able  to  resolve  the  paradox  concerning 
the  shift  of  the  spectrum  with  direction  of  emission. 

Other  significant  contributions  to  the  theory  of  jet  noise 
are  the  discoveries  of  the  correct  nonsingular  convection  factor 
[(1  -  ^  cos  0)2+  a2  Me2  ]“V2  by  Ribner  (Ref.  16)  and  Williams  (P.ef.  17)  and 
of  the  power  laws  {f2  and  f”2)  for  the  total  noise  spectrum  by  Powell  (Sef.l8). 


1.2  Present  Investigation 

The  present  investigation  took  a  more  direct  approach  to 
estimating  the  characteristics  of  the  noise  emission  from  a  unit  volume  of 
turbulence  in  the  mixing  region  of  a  jet  with  very  few  generalizations  and 
assumptions.  Instead  of  using  the  original  Lighthill 's  equation,  the  Proudman 
form  was  used  as  the  basis  of  the  present  work.  This  approach  involves  much 
simpler  experimental  procedure  because  of  the  fact  that  the  noise  intensity 
in  a  certain  direction  is  determined  by  the  velocity  component  at  that  direc¬ 
tion.  For  directions  making  angles  smaller  than  70°  to  the  jet  axis,  the 
velocity  components  can  be  measured  accurately  with  inclined  hot-wires.  Using 
the  original  Lighthill' s  equation  and  the  response  of  inclined  hot-wires, 
Williams  (Ref.  19)  arrived  at  the  same  conclusion. 

Following  Ribner 's  notion  of  'shear'  and  'self  noise,  tvo- 
point  space-time  correlations  of  both  the  turbulent  velocities  and  the  square 
of  these  velocities  were  measured  in  the  mixing  region  of  a  4  inch  low-speed 
model  jet  (about  142  fps)  for  three  directions  9  -•  0°,  45°,  and  6o°,  Pre¬ 
vious  experiments  were  concentrated  on  the  u-components  only.  Using  the  mov¬ 
ing  axis  transformations  of  these  space-time  correlation  functions  to-gether 
with  the  mean  velocity  measurements,  an  attempt  was  made  to  estimate  in  a 
direct  fashion  the  basic  directivity,  the  intensity,  and  also  the  spectrum 
of  both  the  'shear'  and  the  'self'  noise  generated  by  a  unit  volume  of  jet 
turbulence.  With  existing  knowledge  about  the  distribution  of  noise  sources 
in  a  jet,  a  quantitative  comparison  between  the  present  work  and  other 
workers'  acoustic  measurements  of  the  overall  jet  noise  was  possible.  The 
agreement  was  reasonably  good. 


II.  THEORY 


Since  different  ideas  have  been  drawn  from  a  number  of  sources 
to  arrive  at  a  convenient  formalism  for  the  present  experimental  investigation, 
it  is  advisable  to  give  a  brief  coherent  treatment  of  the  theory  here.  Further¬ 
more,  during  the  course  of  this  work,  the  author  came  across  a  simpler  way  of 
expressing  the  spectral  aspect  of  the  acoustic  radiation  problem.  Although  it 
is  mathematically  less  elegant  than  seme  of  the  previous  work,  it  embodies 
most  of  the  main  features  of  this  'aerodynamic  noise'  problem  in  simpler 
fashion.  This  is  given  in  greater  detail  in  Section  2.2. 
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Intensity  Analysis  of  the  Acoustic  Field 


In  this  section,  a  briel  derivation  of  the  acoustic  'far  field' 
pressure  generated  by  a  finite  region  of  turbulent  flow  is  given  using  tensor 
notation  with  the  summation  convention.  Starting  from  the  Navier-Stokes 


equation 


(2.1.1) 


and  the  equation  of  continuity, 


St  +  sl[  (pvi'  =  0  (2.1.2) 

for  an  arbitrary  continuous  fluid  without  sources,  sinks,  or  external  forces, 
Lighthill4  (Ref.  l)  has  obtained  the  following  expression 


where 


&  .  a  2  =  sl¬ 
at2  5xi2 


ss  0  V  V  -  T 

p  i  J  ij 


+  (p  -  an2p)  64 


(2.1.3) 

(2.1.4) 


Now,  if  the  small  effects  of  heat-conduction  and  viscosity  are 
neglected  in  a  medium  without  flow  where  v^,Vj  are  the  velocities  associated 
with  acoustic  waves  only,  the  right-hand  side  of  Eq.  2.1.3  is  negligible  and 
we  have  the  classical  wave  equation  for  a  stationary  acoustic  medium  without 
sound  sources.  On  the  other  hand,  a  medium  with  flow  is  equivalent  to  a 
stationary  medium  containing  a  spatial  distribution  of  virtual  sound  sources 
whose  strength  per  unit  volume  is  given  by  dx^dxj  .  In  other  words 

an  aerodynamic  flow  provides  a  forcing  term  for  generating  an  acoustic  field. 
Because  of  reasons  connected  with  the  double  differentiation,  Lighthill  calls 
the  quadrupole  strength  per  unit  volume. 

Let  a  volume  element  in  the  source  field  be  d^Y  at  point  Y. 
Then  the  solution  for  the  density  field  at  point  X  and  time  t  in  an  unbounded 
medium  is  given  by  the  Kirchhoff  retarded  potential- solution  (Ref.  20) , 


AP  (X,t)  = 


4Tr|x-Y|a02 


SyiSyj 


/  l£-ll  \  , 


(2.1.5) 


where  the  function  5(t')  is  the  Dirac  delta  function.  In  other  words,  the 
quantity  ^T^j/dyjdyj  has  to  be  evaluated  at  the  retarded  time  t  -  |X- Y  |  /aQ. 
In  the  acoustic  'far  field'  where  the  distance  jx|  is  large  compared  with  a 
typical  wave  length  and  the  dimensions  of  the  flow  region,  the  solution  can 
be  rwritten  as  (e.g.,  Refs.  21  and  22). 
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Since  it  is  customary  to  talk  about  acoustic  pressure  in  the 
'far  field',  we  will  rewrite  the  above  equation  in  terms  of  pressure  using 
the  isentropic  relationship, 


Ap(X,t)  - 


V.i 


W;  |x|3 


II 


92 


at‘ 


i' 


t  + 


|X-Y 

TT’ 


-j  dt'  ah 

'  ~  (2.1.7) 


The  mean  square  pressure  is  then  given  by  (Ref.  l), 


I6^a«  |X|‘ 


d3y 


II 


*VZ> 


at. 


0 


|X-Y 


6  {  tg  -  t  +  a 


-) 


dt  d3Z 


(2.1.8) 


where  overbar  denotes  time  average.  For  a  statistically  stationary  turbulent 
motion,  Eq.  2.1.8  can  be  written  as  (Refs.  16,  17,  and  23), 

55  (x) «  rnr  * 


l6w\  |x| 


^  Ti,i©V£’T))<T  -  T*)dTd3Id35 


where 


T  mt2mtl 


T*  = 


„  lx-ll  |x-z| 


(2.1.9) 


o  o 

Equation  2.1.9  can  be  generalized  to  give  the  pressure  autocorrelation  function 
by  introducing  an  arbitrary  time  delay  t',  i.e., 

ap(x)ap(x;^  -  & j]f[ 4?  Tij(^)  vz-t)]  6't  -(■* +  t'” 

n  '  oo 


dTd3^  d3Z  (2.1.10) 


A  considerable  simplification  results  on  expressing  the  tensor  contraction  in 
the  Proudman  form  (Ref.  ll) , 


x.x  T  .  =  |X|2  T 
i  j  ij  xx 


(2.1.11) 


where  is  the  value  taken  on  by  when  the  1-axis  is  aligned  with  the 
vector  X.  Then  we  have, 


Ap(X)Ap(X,T')<* 


TXx(DTxx(Z^) 


5  [t 


d3yd3z 


(t*  +  T  )dT 

(2.1.12) 
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Following  Ref.  l6,  we  introduce  the  following  new  variables  for  experimental 
convenience,  see  Fig.  1, 

H  =  \  (1  +  z) 

=  Y  -  Z  (2.1.13) 

The  corresponding  pressure  autocorrelation  function  reads, 

Ap(X)Ap(X,T  • )  *  -o*  +  -•)]  * 


d3id3ji 


where 


R(H>i>T)  »  +  1/^,t)Txx(I1  -  i/2)  (2.1.14) 

and  the  difference  in  emission  times  given  approximately  by  (Refs.  23  and  24), 


1  •  x 
^!T 


(2.1.15) 


For  reasons  fully  explained  in  his  papers  (Refs.  2  and  25), 
Lighthill  suggests  the  use  of  a  moving  frame  in  analysing  the  turbulent  field 
when  the  eddies  are  moving  at  non-negligible  velocities.  This  is  not  a  re¬ 
quirement  -  as  pointed  out  by  Ribner  in  Ref.  3  -  if  the  differences  in  re¬ 
tarded  time  are  maintained.  However,  for  experimental  convenience  in  neglect¬ 
ing  the  retarded  time  differences,  the  moving  frame  approach  was  chosen  here¬ 
in.  Referring  to  a  co-ordinate  system  moving  with  a  uniform  eddy  convection 
velocity  Uc,  the  pressure  autocorrelation  function  takes  on  the  following 
form  (see  Ref.  17  for  detailed  transformation), 


Ap(X)Ap(X,r • )  * 


l6772a  4  |x  |2  e5 


+  t;  ’] 


drd3i  d3q 
— m  — 


where 


(2.1.16) 


Vi-^c  T 


is  the  equivalent  separation  vector  defined  by  the  transformation  allowing 
the  moving  frame  correlation  function  to  be  defined  in  terms  of  the  fixed 
frame  correlation  function  and  the  time  delay  in  terms  of  _£m,  i.e., 

Rjlb  lm’T  =  R(nd>T) 


-¥•  4.  »  = 

T  ~  T 

'  m  'm 


im-X 


aojX|(1  -  Mccos  9) 


kl  -  Mccos  0, 


8  =  (l  -  M  cosQ  ) 


(2.1.17) 
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ing  relation, 


The  acoustic  intensity  can  then  be  obtained  from  the  follow- 


gg(x) 

Poao 


(2.1.18) 


2.2  Spectral  Analysis  of  the  Acoustic  Field 

Although  a  bit  more  information  regarding  the  power  spectral 
density  can  be  extracted  from  the  four-dimensional  Fourier  analysis  as  used 
by  Kraichman  (Ref.  26),  Mawardi  (Ref.  27),  Lilley  (Ref.  12),  and  Williams 
(Ref.  17),  we  have  found  that  the  essential  features  can  more  easily  be  obtain¬ 
ed  by  the  one-dimensional  Fourier  cosine  transform  technique.  According  to 
the  well  known  Wiener-Khinchin  relations  for  stochastic  processes,  the  power 
spectral  density  of  the  acoustic  pressure  field  is  given  by  the  Fourier 
cosine  transform  of  its  autocorrelation,  i.e., 


(X)Ap(X,Tf )  cosurr'd-r' 


(X)Ap(X,T* )  coscot 'dT' 


(2.2.1) 


Substituting  (2.1.16)  in  (2.2.1), 


'  ITA  -IxlV  /C0S“T'dT,JJi[  i?  >_ 


Changing  the  order  of  integration, 


=  \  - pT  [  d^H  (  l  TT I  Im>T)dl 

l6A.4|X|Vi  J  mJ  St4 


J6  f<T-V)  -  S  ] 


(2.2.2) 


~kZ&J  J  a3-m  J 


l6?ra  |xj  e 


00  00 


Now  defining 


(o,  =  (o  8 
t 


and  expanding  the  cosine  term,  we  have 


cos<o,  (t-t  *)  =  costo,  tcosu),  t  *  +  sinco,  Tsinio,  t  * 
t'm'  t  tm  t  tm 


(2.2.3) 


(2.2.4) 


(2.2.5) 
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Experimental  evidence  (Ref.  28  and  29)  indicates  that  Rm(j[l>.£mT)  is  an  even 
function  in  t,  so  the  sine  term  will  integrate  to  zero. 

Therefore 

°°s  v«*  J  *5  Vt) 


X  COS  U  T  dT 


(2.2.6) 


Since  Rm(il,im,T)  is  even, 

..fen  =  =  0  at  t  =  0 

8t 


If,  in  addition, 


dn*m 

8t  n 


=  0  at  t  =  oo  for  n  -  1,2,3 


then  according  to  the  theory  of  Fourier  transforms  (Ref.  30)  we  have 

1  J*  r 

M-,U)  =  reA^Txl5  fd3]1f  d3lmCOSWt rm*f  \(*VT)cosVdT 


OO  oo 


l67r3a  4  |X|: 


[  d3im  costoer m*J  Rm(jl,lta,T)  cosojGt  dr 


(2.2.7) 


Several  interesting  observations  can  be  made  of  the  above 
equation.  Firstly,  we  can  identify  as  the  frequency  in  the  turbulence;  the 
corresponding  radiated  frequency  ui  is  then  the  Doppler- shifted  frequency 
(i.e.,  w  =  uj^/Q)  .  This  is  a  logical  result  which  one  could  have  obtained  on 
physical  grounds.  Secondly,  Lighthill's  criterion  for  neglecting  retarded 
time  shows  up  automatically  in  the  cosu@  rm*  term.  According  to  Lighthill, 
retarded  time  can  be  neglected  if  ui/a0  is  small  so  that  the  eddy  size  i  is 
small  compared  with  the  wave-length  of  the  sound  it  generates.  If  this  con¬ 
dition  is  met  in  our  case,  then  for  <  £  the  term  cos  can  be  approxi¬ 

mated  as  unity.  Thirdly,  if  retarded  time  is  neglected,  Lighthill's  con¬ 
vection  factor  (l  -  Mc  cos9)"^,  which  accounts  for  the  main  effect  of  ^ 

convection  at  limited  speeds  is  exhibited  as  a  vertical  shift  (l  -  M<>  cos  9) 
plus  a  Doppler  shift  (l  -  Mccos  9)-^  in  the  power  spectrum  and  is  illustrated 
in  Fig.  2.  This  concept  is  a  low- speed  version  of  a  similar  idea  presented  in 
Ref.  3*  Fourthly,  although  the  Light  hill’s  convection  factor  is  not  applic¬ 
able  for  high-speed  convection  because  of  its  singularity  where  1  -  McCos9  =  0, 
the  following  example  will  show  that  Eq.  2.2.7  is  still  valid  for  high  speed 
if  retarded' time  is  not  neglected.  In  fact  this  moving  frame  integral  with 
proper  account  of  retarded  time  possesses  a  zero  that  exactly  cancels  the 
(l  -  MjjCos  9)~5  singularity  and  replaces  it  by  a  nonsingular  convection 
factor  (c.f.  also  Refs.  l6  and  17). 

The  correlation  function  chosen  for  the  present  example  is 
the  equivalent  nonconvective  form  of  the  one  used  by  Ribner  (Refs.  16  and  2k) 
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for  a  homogeneous  and  isotropic  turbulent  flow,  i.e., 


KmtlU.O  -  <*xx>2* 


•  jjl  (tB*  +  «4+  5m3)  +  “f  t2J 


(2.2.8) 


where  L  and  refer  to  a  typical  length  scale  and  frequency  respectively. 
Because  of  the  assumption  of  isotropy,  it  suffices  to  consider  X  to  be  in  the 
1-2  plane.  The  exact  integral  can  then  be  evaluated  directly  with  the  proper 
retarded  time  given  by 


5mi  cos  0  +  5m2  sine 
aQ(l  -  Mc  cos  0) 


(2.2.9) 


Substituting  (2.2.8)  and  (2.2.9)  in  (2.2.7),  we  obtain  for  the  spectrum  of 
noise  from  a  unit  volume  of  turbulence, 


jfeSg)  (*xx)2  L3  a>U 

8%  l67r5/2a04  |x 


Uwf2 


where 


ojf2L2  I 

C  -  [(1  -  M  cos  0)2  +  ]  2 

TTBrf 


(2.2.10) 


(2.2.11) 


The  mean  square  pressure  is  obtained  by  integrating  Eq.  2.2.10  with  respect 
to  frequency, 


d.3Ap2(X)  _  f  a3  (Jj(x,  u)  dw 

J  8  ti 


(2.2.12) 


=  3  L3 

4  V2  a  4  (X|2  C5 

This  result  is  the  same  as  that  obtained  by  Ribner  (c.f.  Eq.  7*9  in  Ref.  3) 
using  a  fixed  frame  formalism  with  a  convective  correlation  function. 

The  above  example  shows  that  when  the  convection  Mach  number 
is  supersonic  the  mean  square  pressure  has  a  finite  peak  at  cos9  =  l/Mc 
instead  of  an  infinite  one,  provided  proper  retarded  time  has  been  considered 
a  fact  first  recognized  by  Ribner  (Ref.  24). 

2*3  'Shear  Noise"  and  'Self  Noise-1 

For  large  Reynolds  number  and  moderate  Mach  number,  T..  can 
be  approximated  by  (Ref.  l),  then 

Txx*  p°vx2  <2*3-D 

where  vx  is  the  instantaneous  velocity  component  along  the  vector  X  from  the 


8 


origin  to  the  observer.  Writing  the  instantaneous  velocity  v  as  the  sum  of  a 
mean  velocity  plus  a  perturbation,  we  have, 

vx  '  Ux  +  ux  <2-3.2) 

and  thus 

R(ji>i>T)  =  P§  [ ux2u^2  +  4UxU^‘uxu^.  +  other  terms]  (2.3*3) 


where  unprimed  quantities  refer  to  the  point  Q  =  +  _|/2  at  time  t  +  t  and 

the  primed  quantities  refer  to  the  point  Q'  =  ^  *  %J 2  at  time  t  (see  Fig.  l). 
The  'other  terms'  are  either  time  independent  and  differentiate  out  or  they 
are  negligible  compared  with  the  first  two  terms.  Ribner  (Ref.  3)  labels  the 
first  term  'self  noise'  and  the  second  term  'shear  noise'. 

Using  a  simplified  model  of  isotropic  turbulence  superposed  on 
a  mean  flow  and  a  joint  Gaussian  probability  density  for  ux  and  u^,  he  obtains 
the  following  relation,* 

4  2 

shear-noise  cos  0  +  cos  0  -  t.\ 

self-noise  2 - 

He  is  also  able  to  obtain  separate  spectra  for  the  'shear'  and  'self'  noise 
peaking  an  octave  apart. 

We  are  going  to  follow  his  notion  of  'shear'  and  'self  noise 
and  try  to  measure  both  the  covariance  of  the  velocity  squared  as  well  as  that 
of  the  velocity  itself.  Thus  we  shall  replace  his  conceptual  model  by  mea¬ 
surements  in  a  real  turbulent  jet. 


Ill .  EXPERIMENTAL  APPARATUS 

3.1  4  Inch  Low-Speed  Model  Jet 

A  detailed  description  of  the  4  inch  model  jet  has  already 
been  given  in  Ref.  28.  Basically,  it  consisted  of  a  screened  diffuser,  a 
settling  chamber  and  a  12:1  contraction  ending  in  a  4  inch  diameter  jet.  It 
was  powered  by  a  one  horsepower  centrifugal  blower  giving  a  velocity  of 
approximately  142  fps.  The  intake  of  the  blower  was  connected  to  a  large 
wooden  chamber  whose  opening  was  covered  with  three  layers  of  ordinary 
furnace  filter  to  provide  a  dust  free  air  flow  (Fig.  3 )• 

3.2  Hot-Wire  Probes  Traversing  Gear 


This  has  also  been  reported  in  Ref.  28.  Essentially,  it  pro¬ 
vided  identical  movement  in  opposite  directions  with  respect  to  a  common  point 
for  the  two  probes.  In  addition,  it  had  three  degrees  of  freedom  for  this 
common  point  (Fig.  3). 


*  Unpublished  work;  this  is  a  correction  of  the  ratio  cos  0  obtained  in 
Refs.  3  and  13. 


9 


J.3  Hot-Wire  Probes  and  Anemometer 


The  hot-wire  probes  were  conventionally  built  from  hypodermic 
.ubing  with  sewing  needles  as  supporting  prongs.  The  wire  itself  was  0.000132 
'nominal  diameter)  tungsten  wire  copper  plated  at  the  ends  to  leave  a  sensing 
ilement  in  the  middle  of  about  l/32  of  an  inch,  A  close-up  picture  of  the 
probes  is  shown  in  Fig.  4. 


The  hot-wire  anemometer  used  was  the  Hubbard  Model  II  HR 
Cype  3A  dual  channel  constant  temperature  anemometer.  It  has  a  built-in 
Linearizing  circuit  and  the  frequency  response  of  the  combination  is  flat  up 
so  ten  kHz  .  As  was  found  in  Ref.  31 >  we  could  obtain  linearity  in  our  speed 
range  only  with  an  overheating  ratio  of  0.3  to  0.4.  The  output  circuit  was 
slightly  modified  to  give  a  smooth  d-c  voltage  reading  for  the  mean  speed. 

?he  d-c  voltages  were  measured  with  two  Simpson  multimeters  which  had  been 
ralibrated  against  a  0.01$  digital  voltmeter. 

J.4  Time-Delay  Correlator 


The  existing  UTIAS  time-delay  correlator  as  reported  in  Ref. 

52  consists  of  essentially  a  two- track  tape  recorder,  an  analog  multiplier,  and 
m  active,  R-C  integrator.  The  adjustable  time  delay  is  obtained  by  a  special 
staggered  play-back  heads  assembly  (Fig.  5)-  In  order  to  meet  the  required 
iccuracy  and  performance  for  the  present  investigation,  the  unit  has  been 
lodified  and  improved.  The  integrator  circuit  was  replaced  by  a  better  per¬ 
formance  (0.2$  accuracy)  active  R-C  integrator  and  the  output  was  read  with  a 
:’our-figures  reading  digital  voltmeter.  An  automatic  controlling  circuit  has 
seen  incorporated  into  the  system.  This  special  circuit  controlled  the  operat- 
.ng  switches  of  the  tape  recorder,  the  starting  and  timing  of  the  integrator 
ind  the  displaying  of  the  result  on  the  digital  voltmeter  in  a  proper  sequence. 
5ee  Fig.  3  for  the  general  set  up.  The  frequency  response  of  the  whole  sys¬ 
tem  was  essentially  flat  from  30  to  10,000  Hz  (Fig.  6).  Its  performance  is 
.llustrated  in  Fig.  7  showing  accurate  correlations  of  sine  waves  of  three 
lifferent  frequencies. 


y.  EXPERIMENTAL  PROCEDURE 

I  The  main  work  of  the  present  investigation  involved  longitudinal 

,wo- point  space-time  correlations  of  not  only  the  turbulent  velocities  but  also 
;he  square  of  these  velocities.  The  measurements  included  the  ordinary 
i-component  and  components  at  45°  and  60°  to  the  jet  axis.  Space  correlations 
>f  these  velocity  components  were  also  obtained  along  the  two  transverse  axes 
,o  give  the  lateral  length  scales.  In  addition,  the  mean  velocity  of  these 
:omponents  were  also  measured. 

-.1  Calibration  of  Hot-Wires 


All  wires  were  calibrated  in  the  laminar  core  of  the  jet  near 
■he  exit  using  a  pitot- static  probe.  The  calibration  was  done  at  the  beginn- 
.ng  of  every  period  of  testing  and  a  two-point  check  was  performed  at  the  end 
•o  ensure  that  there  had  been  no  appreciable  change  in  the  calibration.  Since 
•he  individual  probe  constants  cancel  out  in  forming  the  nondimens ional 
orrelation,  no  attempt  was  made  to  match  the  probe  characteristics.  A  typica] 
alibration  curve  is  shown  in  Fig.  8. 
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4.2  Determination  of  Hot-Wire  Position 

All  our  measurements  were  taken  about  the  location  rjj_/D  =  4 
and  ^2/^  =  0*5  &s  a  center.  The  traversing  mechanism  enabled  the  two  hot-wire 
probes  to  be  moved  at  equal  distances  from  this  point  along  three  orthogonal 
axes.  This  center  point  of  investigation  was  located  by  means  of  two  tele¬ 
scopes:  one  placed  downstream  and  the  other  across  the  jet.  Since  it  was 
fairly  difficult  to  relocate  the  probe  after  a  wire  failure  during  a  run,  that 
particular  test  was  repeated  using  a  new  wire.  However,  such  failure  did  not 
occur  very  often. 

4.3  Determination  of  Hot-Wire  Inclination 

For  measurements  of  the  4-5°  and  60°  components  a  dynamic 
determination  of  the  wire  inclination  was  used  rather  than  a  static  geometrical 
one.  The  wires  were  first  placed  in  the  laminar  core  of  the  jet.  They  v/ere 
then  rotated  to  „ne  proper  inclination  as  indicated  by  the  mean  velocity  mea¬ 
surements  which  were  computed  from  the  usual  sine-law  response  of  inclined 
wires.  (A  discussion  of  the  sensitivity  of  inclined  wires  is  given  in 
Appendix  A.)  The  wires  were  tr.en  transferred  to  the  location  of  investigation. 
Our  special  traversing  mechanism  enabled  such  a  transfer  without  changing  the 
inclination  of  the  wires  to. the  jet  axis. 

4.4  Method  of  Measurements 

For  each  separation  of  the  two  hot-wire  probes  in  the  flow, 
the  electrical  signals  that  were  related  to  the  turbulent  velocity  components 
were  recorded  simultaneously  on  the  two  tracks  of  the  tape  recorder.  The  d-c 
readings  that  were  proportional  to  the  mean  velocities  were  recorded  separately. 
When  correlating  the  square  of  turbulent  velocities  the  signals  were  first 
squared  by  means  of  a  dual-channel  Fhilbrick  Analog  Multiplier  before  they 
were  taped.  At  the  initial  position  of  zero  separation  the  true  rms  values 
of  the  a-c  signals  were  measured  with  a  Flow  Corporation  Random  Signal  Meter. 

The  product  of  the  two  wires'  readings  served  to  give  the  mean  square  values 
of  the  turbulent  velocities  at  the  center  point.  After  readings  at  maximum 
separation  were  taken,  the  two  probes  were  returned  to  their  initial  positions 
and  a  check  was  performed  on  both  the  a-c  and  d-c  readings.  The  deviation  was 
usually  less  than  Elock  diagrams  of  the  recording  and  data  processing 

system  are  shown  in  Figs.  9  and  10. 

In  order  to  detect  any  changes  in  sensitivity  of  the  record¬ 
ing  tape  and  possible  electronic  drift,  a  fixed  intensity  1,000  Hz  sine  wave 
was  recorded  at  fixed  intervals  during  the  run  to  serve  as  reference  points. 

Such  possible  changes  were  later  corrected  by  assuming  a  linear  variation 
between  the  reference  points  in  the  data  processing. 

When  the  two  recorded  signals  were  played  back,  a  finite  time 
delay  t  was  inserted  between  them  by  means  of  the  adjustable  staggered  playback 
heads.  Then  the  signals  were  amplified,  multiplied,  and  averaged  by  means  of 
an  active  R-C  integrator  whose  integration  time  is  10  seconds.  The  results 
were  normalized  with  respect  to  the  mean  square  value  at  the  center  point. 

This  process  was  repeated  for  a  sequence  of  values  of  t  . 
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V.  EXPERIMENTAL  RESULTS 


The  main  objective  of  the  present  investigation  has  been  to 
estimate  the  intensity,  the  spectrum,  and  the  basic  directivity  of  the  noise 
generated  by  a  unit  volume  of  turbulence  fr^ji  the  two-point  space-time  corre¬ 
lation  measurements  of  turbulent  velocities.  However,  we  feel  it  is  worth 
reporting  here  first  some  of  the  characteristics  of  jet  turbulence,  particularly 
the  45°  and  60°  components  which  have  not  been  measured  before.  Whenever 
possible,  the  present  results  are  compared  with  those  of  other  workers. 


5.1  Mean  Velocity  Measurements 

The  variation  of  the  mean  velocity  across  the  jet  has  pre¬ 
viously  been  measured  with  a  single  hot-wire  by  the  present  author  and  report¬ 
ed  in  Ref.  33*  A  similarity  parameter  has  been  rediscovered  which  enabled  the 
mean  velocity  profiles  of  jets  up  to  low  supersonic  speed  to  be  plotted  on  one 
single  universal  curve. 

r 

5.2  Turbulence  Intensity 

A  complete  survey  of  the  turbulent  intensities  across  the  jet 
has  not  been  carried  out  herein.  But  parts  of  the  turbulent  intensity  pro¬ 
files  around  t^/D  =  0.5  were  obtained  from  recorded  data  of  two-point  space- 

time  correlations.  They  are  shown  in  Fig.  11(a).  Note  that  the  intensity  of 
the  45°  component  was  higher  than  that  of  the  ordinary  u-component,  and  that 
of  the  60°  component  was  lower.  Other  workers  (Ref.  9)  have  shown  that  the 

intensity  of  the  90°  (or  v)  component  was  lower  than  that  of  the  u-component, 

which  seems  consistent.  Figure  11(b)  shows  a  favourable  comparison  between 
the  present  result  and  that  of  Laurence  (Ref.  7)  for  the  nondimensional 
quantity  /u^/(/u^)max  at  the  same  downstream  location. 


5.3  Space  Correlations  and  Length  Scales 


The  signals  from  two  hot-wire  probes  separated  by  a  distance 
?  can  be  correlated  to  define  a  region  over  which  the  velocities  at  the  two 
points  are  strongly  related.  These  results  can  then  be  used  to  define  an 
integral  length  scale  or  typical  'eddy'  size.  For  conveniences  in  later 
computational  work,  the  space  correlation  coefficient  is  defined  as  follows, 


(R  (1  ,i,e)  = 


u9  (H  +  1/2)  ue  in  -  i/2) 

tn) 


(5.3.1) 


where  uq  is  the  turbulent  velocity  component  making  an  angle  0  with  the  jet 
axis *and  the  overbar  denotes  an  average  over  a  period  of  time  sufficiently 
long  to  obtain  a  stationary  value. 


The  integral  length  scale  is  defined  as 
•  00 

Lug  (a)  =  J  (Rind**)  di 

or  alternatively  as  ’ 

4^(2)  =  f~iRin>l>s)d 

0 


(5.3.2) 

(5.3.3) 
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depending  on  experimental  convenience  (see  later),  where  is  the  maximum 
separation  at  which  $  (^,.1,9)  goes  to  zero.  Similar  definitions  are  used  for 
the  correlations  of  the  square  of  the  velocities. 

The  nine  quantities  <R  (3,^, 6)  are  shown  in  Fig.  12  for 
correlations  of  both  the  turbulent  velocities  and  the  square  of  these  velocit¬ 
ies.  The  general  forms  of  the  u-component  correlations  in  the  three  ortho¬ 
gonal  axes  were  in  close  agreement  with  those  of  Ref.  9*  The  decrease  in 
magnitude  of  the  negative  correlation  with  increasing  0in  the  circumferential 
traverse  was  also  compatible  with  results  of  Ref.  9  which  showed  a  non-nega¬ 
tive  correlation  for  the  component  at  8  =  90°. 

In  computing  the  integral  length  scales,  Eq.  5*3*2  was  used 
for  the  radial  and  circumferential  cases  whereas  Eq.  5*3*3  was  used  for  the 
axial  case.  In  this  connection,  the  negative  correlation  in  the  axial  tra¬ 
verse  was  most  likely  a  spurious  effect  due  to  inadequate  low-frequency  cut¬ 
off  of  the  tape  recorder  in  the  time-delay  correlator.  A  detailed  discussion 
of  the  effect  of  low-frequency  cut-off  is  given  in  Appendix  B.  Estimates 
based  on  the  finding  in  Appendix  B  indicated  that  our  transverse  length  scales 
might  be  lower  by  about  10$  and  that  the  longitudinal  scales  might  be  off  by 
as  much  as  20$.  The  different  length  scales  nondime:  Lonalized  by  the  jet 
diameter  are  given  below. 

for  correlation  of  Uq 


(«!»  0,  0)* 

(0,  S2,  0) 

(0,  0,  s3) 

V* 

0.191 

0.138 

0.023 

kul^o/0 

0.168 

0.187 

0.025 

Nso^13 

0.144 

0.223 

0.048 

2 

for  correlation  of  u_ 

iiv  0,  0)* 

(0,  s2,  0) 

(0,  0,  S3) 

0.102 

0.057 

0.050 

I*u.2i^o/D 

0.080 

0.059 

0.041 

0.070 

0.061 

0.044 

The  small  length  scale  associated  with  the  correlation  of  uq 
along  the  axis  ('third  column  of  first  table)  arises  from  the  large 
negative  loops  in  Fig. 12.  This  does  not  imply  that  the  eddies  are  physically 
very  short  in  the  direction. 
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5.4  Two-Point  Space-Time  Correlation 


If  the  signal  from  the  upstream  probe  is  delayed  a  time  t 
before  being  correlated  with  the  other,  the  two-point  space- time  (or  cross) 
correlation  will  be  obtained 


_  ,  „  »,(a  ♦  WF)  .(a  -  W> 

-4(a) 


(5.4.1) 


A  similar  definition  is  used  for  correlations  of  the  square  of  velocities. 
Results  for  the  three  components  in  the  axial  traverse  are  shown  in  Figs. 

13  -  18.  The  convective  nature  of  the  turbulent  field  is  revealed  by  the 
’yawing'  of  the  correlation  'ridge'  of  the  typical  three  dimensional  plot  in 
Fig.  19  (e.g.,  Ref.  16). 


The  correlation  coefficients  for  different  probe  separations 
did  not  rise  to  unity  with  optimum  time  delay  as  shown  in  previous  figures. 
This  fact  showed  that  the  turbulence  'pattern'  was  not  a  'frozen  pattern'. 
According  to  Ref.  29  the  envelope  of  these  cross-correlation  curves  is  the 
autocorrelation  of  the  fluctuation  which  would  be  seen  by  an  observer  moving 
with  the  turbulence.  The  Fourier  cosine  transform  of  this  evenlope  will  give 
the  spectral  density  of  the  turbulence  relative  to  axes  moving  with  the  flow. 
Figure  20  shows  a  comparison  between  the  fixed  frame  and  the  moving  frame 
spectra  of  the  axial  turbulent  velocity.  The  fixed  frame  spectrum  was 
obtained  as  usual  with  a  Murhead  wave  analyser.  The  moving  frame  spectrum 
was  obtained  from  the  moving  frame  autocorrelation  using  Filon's  numerical 
method  (Ref.  34).  The  computation  was  done  with  the  IBM  7094  computer  at 
the  Institite  of  Computer  Science  in  the  University  of  Toronto.  It  is  seen 
that  the  fixed  frame  spectrum  contains  more  of  the  high  frequencies  due  to 
the  convection  of  a  spatial  pattern  passing  the  fixed  hot-wire  probe. 


Richards  and  Williams  (Ref.  35)  considered  this  moving  frame 
turbulence  spectrum  and  speculated  that  the  radiated  sound  spectrum  could 
roughly  be  deduced  from  it  by  the  addition  of  twelve  db  per  octave.  This  is 
equivalent  to  multiplying  the  moving  fVame  spectrum  by  w  .  However,  the 
very  strong  amplification  of  the  high  frequencies  makes  accurate  estimation 
very  unlikely.  Also  the  assumption  that  the  simple  u-velocity  spectrum 
will  represent  the  spectrum  of  the  more  complex  T^  is  a  great  over¬ 
simplification.  Therefore,  there  is  no  escape  from  using  a  more  elaborate 
approach  as  given  in  Section  VII  here. 

5.5  Convection  Velocity 

For  a  frozen  pattern  as  found  in  turbulence  behind  grids  the 
turbulent  convection  velocity  is  equal  to  the  mean  flow  velocity.  In  the 
case  of  turbulent  shear  flow  the  shear  stresses  distort  the  convected  pattern 
as  it  travels  downstream.  As  a  result  there  is  no  unique  convection  velocity. 
However,  from  the  view  point  of  aerodynamic  noise  generation,  a  good  criterion 
for  the  definition  of  convection  velocity  is  that  the  time  scale  of  the 
fluctuating  turbulent  stress  should  be  a  maximum  in  a  reference  frame  moving 
with  this  velocity.  Thus  a  meaningful  definition  of  the  convection  velocity 
is  the  ratio  ?i/t  at  the  point  where  (R  (ij, I-^jt)  =  0  .  As  has  been 

pointed  out  by  Williams  (Ref.  29)  and  Fisher  and  Davies  (Ref.  36).  This 
corresponds  to  the  convection  velocity  defined  by  the  time  delay  for  which 


14 


the  curve  of  (ft  (t^,  t)  versus  t  coincides  with  the  envelope  of  all  such  curves 
From  Figs.  13  -  18,  the  following  values  for  convection  velocity  were  obtained; 


Cases: 


uo 

u45° 

CM 

O 

3 

O 

O 

3° 

u45o 

u2 

u6o° 

0.58 

0.61 

0.6l  0.62 

O.63 

0.66 

Within  experimental  error,  it  seemed  that  the  different  velocity  components 
travelled  with  the  same  convection  velocity.  This  was  in  agreement  with  re¬ 
sults  of  Ref.  9  for  the  axial  and  radial  components.  Also  the  approximate 
magnitude  was  in  agreement  with  results  of  Refs.  8  and  9  at  the  same  cross 
stream  location. 


VI.  ESTIMATION  OF  NOISE  INTENSITY  AND  BASIC  DIRECTIVITY 

In  our  case  of  a  low  velocity  jet,  the  differences  in  retarded 
time  within  a  correlation  volume  could  be  neglected  because  the  'eddy'  size 
was  small  compared  with  the  typical  wavelength  of  sound  it  generated.  A 
posteriori  check  showed  that  the  ratio  was  about  0.01.  Neglecting  retarded 
time,  we  can  interchange  the  integration  and  differentiation  processes,  and 
Eq.  2.1.16  then  becomes 


Ap(x)Ap(x,T‘)* 


l67T2a^  |x|2c5 


! 


6(t~ 


T  ’)dT 

m  ' 


(6.1) 


Note  that  the  Lighthill’s  convection  factor  6”^  has  been  replaced  by  the  non¬ 
singular  factor  C-5.  This  may  be  justified  as  follows.  We  have  seen  in 
material  leading  to  Eq.  2.2.12  that  when  retarded  time  is  included  the  final 
convection  factor  is  of  the  cited  form  C"  .  Therefore  the  replacement  of  0~5 
by  C“5  here  should  constitute  a  first  approximation  to  a  correction  for  the 
neglect  of  retarded  time.  The  intensity  of  noise  generated  from  unit  volume 
of  turbulence  follows  from  Eq.  2.1.18  as 


d3l(X) 


_ 1 _ 

l6772PGao5  l£l2c5 


(6.2) 


In  the  present  co-ordinate  system  there  was  convection  only 
in  the  axial  (or  l)  direction,  therefore,  from  Eq.  2.1.17, 


»  ^mp*  ^mg’  ~  ^(jlj  +  6cT,f2>£3>T)  (6.3) 

Note  that  even  for  the  simple  case  of  estimating  R  (tj, Ojsn^jO,-!-)  one  of  the 

probes  has  to  be  displaced  not  only  in  the  2-direction  but  also  in  the 
1-direction  by  an  amount  of  Uct.  The  fact  that  the  convection  velocity 
varies  across  the  jet  (Ref.  8  and  9)  adds  more  complication  to  the  investi¬ 
gation.  Thus  to  hold  the  amount  of  labour  within  reasonable  bounds  we  were 
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forced  to  use  the  simplifying  assumption  that  the  four-dimensional  space- 
time  correlation  function  was  separable  in  space  in  the  two  transverse  axes 
while  keeping  the  nonseparable  nature  in  space  and  time  in  the  axial  direction, 

i.e^ 

R(n»*mi  +  Uct»I2»53>t^  a  Rito,5iai  *  UcT»T)  R2^»52^  r3^»53)  (6**0 

Note  that  the  present  assumption  is  less  restrictive  than  the  usual  one  of  a 
completely  separable  function  in  space  and  time  (Refs.  3>  12,  and  17).  With 
this  assumption,  it  can  be  shown  that 

HQ.L*  V-S2>  V> - 1 - ° - 2 - - - L~ 

13  [S(J!,0,0,0,0)]2  (6_5) 


The  intensity  per  unit  volume  becomes, 


a3i(x)  ^  R(j,o, o,o,o)  r00  r(h,o,o,s3,o)  r' 

d3H  27^0  O  |X|2C^  R(jl»0,0,0,0)  ^  Jq 


^(,3»0>  52>0>o)  ^ 
R(t!,0, 0,0,0)  2 


r  e(T)dT  r  r  Rxn>gmi  +  uct>  q>q>t) 

J  L  o  R(n,o ,0,0,0) 


3ml  I 

(6.6) 


Using  results  of  Section  2.3,  we  have  for  the  ’ shear ’noise: 


2p0ux2(a)u?£a)  r"  ,,  r  du  "  shB(ji,5mi+ucT,o,o,T) 

L;  ^ixi2 <*  -  z&l6(rn^l 


where 


Ux^x’  uxux  (jl>0>^5 ?3>0) 


ux2(il) 


l®(h>5]_>0,0,t)  =  UXUX  uxux(jj,  Ii,0,0,t) 


(6.7) 


and  for  the  'self*  noise: 


°  3  self 


-5  1,1s  /B(-r)dxL|!_  f SeB  tidS]*  d|„ 


where 


a 

-X 


ux2ux2  (h>0>0553j0) 

—  =*TT  —  d^3  et< 

ux  'i5  _ 

;R(H>il>0>0>T)  =  UX2U'2  (j}>  £^>0,0,t) 


(6.8) 
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Knowing  the  convection  velocity,  R(Tj,|m  +  U  t,0,Qt)  can  be  obtained  from 
R(n.t1.0,0,r). 


In  obtaining  for  the  'shear'  noise,  the  mean 

velocities  were  first  corrected  for  the  effect  of  turbulence  using  the  follow¬ 
ing  equation  (Ref.  37) » 


U 


corr. 


U 

exp. 


exp. 


(6.9) 


Figures  21  -  26  show  plots  of  UxtP 

Ux2U)ux<n) 

for  the  three  components  in  the  two  transverse  axes  together  with  the  re¬ 
spective  areas  under  these  curves.  This  quantity  in  the  axial  direction  will 
be  the  same  as  the  correlation  function  defined  in  Sec.  5*^  because  of  the 
uniformity  of  the  mean  velocities  within  our  correlation  separations.  Figures 
27  and  28  show  the  typical  forms  ot  R(jj,!mi+  Uct,0,0,t)  as  derived  from  the 
fixed  frame  correlation  for  both  tne  'shear'  and  'self 'noise.  As  has  been 
pointed  out  earlier,  the  negative  correlation  at  large  separation  was  most 
likely  due  to  the  inadequate  low  frequency  response  of  the  instrument.  Re¬ 
sults  from  Appendix  B  suggest  that  to  a  first  approximation,  the  correction 
term  due  to  low-frequency  cut-off  when  transformed  to  the  moving  frame  is 
independent  of  time  delay.  This  conclusion  is  manifested  by  the  collapsing 
of  the  families  of  curves  onto  a  single  curve  at  large  separation  shown  in 
Figs.  27  and  28.  The  integrated  areas  of  R(jj,lm1+  Uct,0,0,t)  versus  £mi  as 
a  function  of  t  are  shown  in  Figs.  29(a)  and  30(a).  The  problem  now  is  to 
determine  the  fourth  derivative  of  these  final  curves  with  respect  to  t. 


6.1  Extraction  of  Higher  Order  Derivative 


Since  the  whole  problem  hinges  on  the  extraction  of  the  fourth 
derivative  of  some  experimental  curves,  a  sophisticated  method  is  required. 

Iu  Ref.  38,  several  methods  have  been  tried  to  obtain  the  second  derivative 
of  experimental  data.  It  was  concluded  there  that  the  traditional  technique 
of  fitting  (by  least  square)  a  high  degree  polynomial  to  the  data  was  still 
the  best  method.  For  the  present  purpose  a  more  elaborate  IBM  709*+  computer 
program  was  written  whereby  the  data  points  could  be  fitted  by  a  polynomial 
of  any  degree.  Only  even  power  terms  in  the  polynomial  were  used  because  we 
were  dealing  theoretically  with  an  even  function. 

It  was  impossible  to  estimate  in  any  simple  ways  the  probable 
error  in  the  final  experimental  data  points,  so  the  criterion  of  best  fit  could 
only  be  furnished  by  the  sum  of  the  standard  derivations  (Bef.  39)*  The  num¬ 
ber  of  terms  in  the  polynomial  were  increased  one  at  a  time  until  the  sum  of 
the  standard  deviations  was  approaching  a  minimum.  It  turned  out  that  the 
fourth  derivative  at  the  origin  was  also  approaching  approximately  an  asymp¬ 
totic  value. 


This  approach  was  first  tested  out  with  a  Gaussian  function. 
This  function  was  sampled  at  approximately  the  same  intervals  as  our  experi¬ 
mental  data.  In  order  to  introduce  some  arbitrary  error  the  magnitudes  of 
these  sampled  points  were  rounded  off  at  two  decimal  places.  Figure  31  shows 
that  the  comparison  between  the  theoretical  fourth  derivative  and  that  obtained 
from  the  best  polynomial  fit  is  very  good  up  to  a  fair  range  of  t. 
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Since  we  do  not  have  the  complete  four-*  derivative  curve 
i'rom  the  polynomial  approximation,  it  is  .ecessary  to  investigate  the  range 
of  the  fourth  derivative  curve  required  to  give  a  proper  Fourier  cosine 
transform  for  estimating  the  spectrum.  It  turned  out  that  when  the  fourth 
derivative  curve  resembled  a  damped  oscillation,  the  first  half  wave-length 
would  suffice  to  give  a  fairly  good  estimate  of  both  the  peak  and  the  general 
shape  of  its  Fourier  cosine  transform.  This  is  illustrated  in  Fig.  32,  still 
using  the  Gaussian  function.  The  fourth  derivative  of  the  approximating 
polynomial  was  fitted  with  an  analytical  function  ^e“Bir2C0, 

whose  Fourier  cosine  transform  was  then  compared  with  the  theoretical  result. 

When  this  technique  was  applied  to  actual  experimental  results, 
two  other  types  of  function  were  tried  in  addition  to  the  straight  polynomial. 
They  were 

-(V  a2T2  + ...  +  Vn> 

(1)  e 


_ 1 _ _ 

p 

aQ  +  a2T  +  . . .  +  anTn 

These  were  obtained  by  fitting  an  even  degree  polynomial,  still  in  the  least 
square  sense,  to  the  logarithm  and  the  reciprocal  of  the  experimental  data 
respectively.  A  comparison  of  the  typical  results  for  the  three  types  of 
fitting  is  given  in  Figs.  33  and  34.  Note  the  consistant  behaviour  of  the 
fourth  derivative  at  the  origin  as  a  function  of  the  degree  in  the  fitting 
polynomial  for  the  three  cases.  It  must  be  pointed  out  here  that  even  if 
the  magnitude  of  the  fourth  derivative  at  the  origin  changes  by  a  factor  of 
two,  the  resulting  noise  intensity  will  change  by  only  three  db.  This  is 
not  considered  to  be  a  serious  error  in  acoustic  measurements.  The  expon¬ 
ential  function  with  a  polynomial  exponent  was  chosen  for  consistency  because 
it  gave  the  damped  oscillatory  shape  of  the  fourth  derivative  for  all  cases 
considered  (Figs.  29(b),  and  30(b)).  The  straight  polynomial  fit  was  re¬ 
jected  because  the  fpurth  derivative  obtained  from  it  resembled  a  pure  cosine 
function.  The  Fourier  cosine  transform  of  it  would  be  a  delta  function  im¬ 
plying  that  a  pure  tone  was  generated  by  the  unit  volume  of  turbulence  con¬ 
sidered.  Physically  this  is  less  plausible  than  the  peaky  spectra  obtained 
from  the  other  two  functions. 

Although  the  least  square  curve  fitting  method  seems  to  be 
the  most  promising  and  consistent  way  of  dealing  with  our  data,  it  should  be 
pointed  out  that  the  results  so  obtained  are  quite  sensitive  to  the  range  of 
data  available.  When  we  first  applied  this  method  to  our  data,  we  found 
that  at  zero  degree  to  the  jet  axis,  the  'shear'  noise  intensity  was  about 
twenty-two  times  higher  than  that  of  the  'self'  noise  and  that  the  peak  fre¬ 
quency  of  the  'self'  noise  was  only  1.3  times  that  of  the  'shear'  noise. 

These  results  are  incompatible  with  measured  jet  noise. 

Upon  carefully  re-examining  the  data,  we  observed  that  be¬ 
cause  of  the  limited  abscissa  scale  we  had  only  about  25$  variation  in  magni¬ 
tude  of  our  correlation  for  the  'shear'  noise  cases  as  compared  to  the  45$ 
variation  for  the  'self'  noise  cases  (Figs.  29(a),  and  30(a))*  In  order  to 
obtain  a  fair  comparison  between  the  two  types  of  noise,  the  'shear'  noise 
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data  should  be  extrapolated  to  give  approximately  the  same  percentage  variation 
in  magnitude  as  those  of  the  'self'  noise.  Exponential  functions  of  the  type 

e"(aT)  were  used  in  the  extrapolation  by  replotting  the  logarithm  of  the 
tail  end  data  of  Fig.  29(a)  on  a  log  log  scale  (Fig.  35).  Extrapolated  values 
of  the  three  'shear1  noise  cases  were  shown  as  black  dots  in  Fig.  29(a). 

From  results  of  Section  VI,  the  length  scales  obtained  in 
Sec.  5 * 3 a  and  the  fourth  derivative  of  the  moving-frame  longitudinal  space- 
time  correlation  at  zero  time  delay  obtained  here,  an  estimate  of  the  noise 
intensity  generated  by  a  unit  volume  of  turbulence  could  be  obtained  through 
Eqs.  6.6  and  6.7.  A  listing  of  the  coefficients  in  the  polynomial  exponent  of 
the  exponential  functions  used  in  the  least  square  curve  fitting  is  shown  in 
Appendix  C.  The  corresponding  curves  are  shown  in  Figs.  29(a)  and  30(a). 

The  final  results  are  given  in  the  following  table  in  the  non- 
dimensional  form. 

d3  /l(x)  #a05  j|foV\ 

a^T)  PqU3  'shear 


a3T) 


C- 


(X)7T2a05  |X|2C5U 


poU 


self 


0°  case  11.48  x  10"5  4.36  x  10“5 

45°  9-36  x  10-5  I.63  x  10" 5 

60°  6.63  x  10-5  3.77  x  10-5 

Note  that  the  results  given  here  are  for  noise  generated  from  a  unit  volume 

of  turbulence  to  an  observer  in  the  far  field  located  in  the  plane  of  the 
jet  axis  and  the  center  of  our  unit  volume.  Therefore  one  cannot  really  in¬ 
fer  the  emission  in  direction  0  from  a  slice  of  jet  without  results  from 
different  unit  volumes  around  the  jet  at  the  same  axial  and  radial  location. 
However  the  differences  are  not  expected  to  be  large.  A  more  detailed  dis¬ 
cussion  of  this  aspect  is  given  in  Appendix  D. 

One  way  to  check  the  consistency  of  our  approach  in  extract¬ 
ing  the  fourth  derivative  is  by  means  of  the  peak  frequencies  obtained  in  a 
later  section.  Since  the  radiated  intensity  is  proportional  to  (fpeak)^ 
approximately,  other  things  being  equal,  the  ratio  of  the  fourth  derivative 
at  zero  time  delay  for  the  'self'  and  'shear'  noise  normalized  correlation 
functions  should  roughly  equal  to  the  ratio  of  their  peak  frequencies  raised 
to  the  fourth  power.  The  results  showed  that  they  were  approximately  equal, 
being  56  and  42  respectively. 

Within  the  limitation  of  our  simplifying  assumptions,  the 
present  investigation  indicated  that  at  zero  degrees  to  the  jet  axis  the 
'shear'  noise  was  about  2.6  times  higher  than  the  'self'  noise  generated  by 
a  unit  volume  of  turbulence  in  the  mixing  region.  Had  we  followed  other 
investigators  (Refs.  12,  13,  and  17)  in  assuming  complete  separability  for 
the  four-dimensional  space-time  correlation  fmiction  and  approximating  the 
moving-frame  autocorrelation  function  by  e"wf’lT'  we  could  have  obtained 
approximately  equal  magnitude  for  the  'shear'  and' self’  noise  as  suggested  in 
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Ref.  13.  However,  a  ratio  of  1.7  was  obtained  if  a  Gaussian  function  was 
used  instead  of  the  exponential  function.  In  view  of  the  compatibility  test 
of  the  last  paragraph,  our  present  results  might  not  be  too  far  from  reality. 

The  basic  directivity  for  the  ’shear'  noise,  as  inferred  from 
the  present  experiment  in  three  directions,  varied  about  like  |cos0|  (|cos0|°*7 
to  be  more  exact)  over  the  range  from  0°  to  60°  from  the  jet  axis.  Since 
UxUj.  =  0  at  90°  it  seems  reasonable  to  extrapolate  the  |cos0 (variation  up  to 
90°.  This  directivity  resembled  better  Ribner's  (cos^0+cos20 )/2  pattern* 
than  the  four-leaf  clover  sin20cos20  suggested  by  Lighthill  (Ref.  2)  and 
Lilley  (Ref.  12).  For  the  case  of  ’self’  noise,  a  minimum  value  was  obtained 
at  45°.  This  was  considered  as  more  likely  an  experimental  error  than  a  true 
physical  situation.  Results  for  0°  and  60°  suggested  on  the  other  hand  that 
the  'self'  noise  had  a  more  or  less  omnidirectional  pattern.  This  was  in 
agreement  with  Refs.  3  and  13. 

The  overall  basic  directivity  would  roughly  be  given  by 
(l  +  2.6|cos0|).  This  result  is  plotted  in  Fig.  36  where  the  convection 
factor  C“5  has  been  incorporated.  Note  that  the  basic  directivity  pattern 
shows  no  valley  along  the  jet  axis.  The  observed  valley  or  dimple  in  all 
acoustic  measurements  &f  jet  noise  is,  in  fact,  due  to  refraction.  This  is 
well  substantiated  by  recent  work  done  at  ITTIAS  (Ref.  40) .  The  effect  of 
refraction  qualitatively  sketched  in  as  the  dashed  curve  in  Fig.  36. 


VII.  ESTIMATION  OF  NOISE  SPECTRUM 


Following  the  assumptions  made  in  Sec.  VI,  Eq.  2.2.6  for  the 
spectrum  function  can  be  written  as  in  the  following, 


for  the  'shear'  noise: 
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and  for  the  'self'  noise: 
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'KXji’J/o  £  refers  to  the  noise  spectrum  from  unit  volume  of 
The  problem  here  is  to  determine  the  Fourier  cosine  transform  of 
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pattern  obtained  in  Refs.  3  and  13 . 
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the  fourth  derivative  curves  obtained  in  the  last  section.  Following  the 
example  with  the  Gaussian  function,  an  analytical  function  (A  sech  Btcoso)0t) 
was  first  fitted  to  the  fourth  derivative  curve  and  an  analytical  spectrum 
function  could  then  be  obtained.  In  addition,  a  numerical  Fourier  cosine 
transform  was  performed  on  the  fourth  derivative  curve  up  to  the  third  zero 
crossing,  assuming  zero  value  thereafter.  Filon's  method  was  used.  Figure 
37  shows  the  approximating  curves  and  Fig.  38  shows  the  corresponding  spectrum 
functions  for  the  axial  direction  case.  The  occurrence  of  negative  values  in 
the  spectrum  function  in  the  second  method  was  due  to  truncation  of  the  fourth 
derivative  curve  at  a  finite  r  as  illustrated  in  Appendix  E.  Thus  it  was 
felt  that  the  spectrum  function  given  by  the  analytical  case  was  the  most 
probable  one,  with  the  uncertainty  given  by  the  shaded  area  in  Fig.  38.  The 
results  indicated  that  the  ' self*  noise  had  a  peak  frequency  which  was  about 
2.5  times  that  of  the  ’shear*  noise.  This  was  even  higher  than  the  factor  2 
inferred  by  Ribner  (Ref.  13)  from  the  squaring  of  the  turbulent  velocities 
for  the  'self'  noise.  However,  the  general  shape  of  the  spectrum  function 
was  in  agreement  with  his  function 

-2  — p 

4 

uk  e 


where  cjp  is  the  peak  frequency. 

The  spectra  of  the  noise  intensity  at  45°  and  6o°  to  the  jet 
axis  would  have  very  nearly  the  same  peak  and  shape  as  for  the  0°  case 
judging  from  the  fourth  derivative  curves.  These  will  not  be  presented. 


VIII.  COMPARISON  WITH  ACOUSTIC  MEASUREMENTS 


Before  the  quantitative  comparisons  are  presented,  one  should 
be  aware  of  the  limitations  and  extra  assumptions  involved.  Firstly,  we 
have  to  infer  emission  for  a  slice  of  jet  from  volume  element  measurements 
taken  in  a  plane  containing  the  observer  and  the  jet  axis.  Secondly,  we  have 
to  rely  on  existing  knowledge  about  the  distribution  of  noise  sources  in  a 
jet  based  on  theoretical  models.  Thirdly,  extrapolation  of  higher  speed  jet 
data  was  used  because  of  the  unavailability  of  acoustic  measurements  at  our 
low  speed  range. 


8.1  Dominant  Frequency 

There  have  been  conflicting  reports  (Refs.  6,  4l,  42,  43,  44, 
and  45)  as  to  the  dependence  of  the  peak  frequency  in  the  noise  spectrum  on  j el 
velocity.  Lee's  results  (Ref.  43)  were  chosen  for  comparison  because  they 
were  in  the  relatively  low  speed  range.  Figure  39  is  a  replot  of  his  results 
taken  at  25°  from  the  jet  axis  and  they  suggested  that  the  peak  frequency 
varied  like  U^*58.  The  variation  would  be  even  smaller  if  Doppler  shift  had 
been  accounted  for.  Using  the  results  (Refs.  46  and  47)  that  the  peak  fre¬ 
quency  was  inversely  proportional  to  the  jet  diameter,  a  peak  frequency  of 
about  330  Hz  was  estimated  for  the  present  4  inch  jet  from  Lee's  results. 

Now  the  noise  spectrum  for  the  whole  jet  can  be  considered  as 
being  made  up  of  spectra  from  different  slices  of  the  jet.  In  the  mixing 
region,  the  spectrum  of  noise  from  a  slice  of  jet  will  have  approximately 
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the  same  peak  as  that  from  a  unit  volume  at  the  region  of  maximum  shear  be¬ 
cause  the  chief  noise-emitting  eddies  are  confined  in  a  small  volume  at  this 
region  (Ref.  12).  Results  of  Ref.  15  suggested  that  the  peak  frequency  of  the 
over-all  noise  spectrum  was  generated  by  a  slice  located  at  about  5  diameters 
downstream.  Similarity  consideration  (Ref.  18)  implied  that  frequency  was 
inversely  proportional  to  axial  distance  in  the  mixing  region.  From  all  these 
considerations,  the  peak  frequency  of  the  noise  generated  by  a  unit  volume  of 
turbulence  at  the  position  of  the  present  investigation  was  estimated  to  be 
about  410  Hz.  The  present  experimental  result  gave  a  value  oi  approximately 
340  Hz  which  was  1.2  times  lower.  This  is  the  value  of  the  'shear'  noise 
spectrum  at  zero  degree  since  the  'shear'  noise  spectrum  dominates  at  that 
direction.  In  view  of  the  discrepancies  of  different  data  and  all  the  uncer¬ 
tainties  in  estimating'  the  slice  peak  frequency,  the  comparison  is  considered 
quite  satisfactory. 

8.2  Reverse  Doppler  Shift  Paradox 

Assuming  that  the  relative  separation  between  the  ’shear’  and 
'self'  noise  spectra  was  the  same  for  different  slices  of  the  jet  in  forming 
the  total  spectra,  our  results  would  predict  qualitatively  one  observed 
'paradox'  in  acoustic  measurements  of  jet  noise.  References  6,  43,  and  44 
reported  that  the  noise  spectrum  at  small  angle  to  the  flow  peaked  at  lower 
frequency  than  that  at  right  angles.  This  is  just  the  reverse  of  the  expecta¬ 
tion  from  a  single  spectrum  subjected  to  Doppler  shift  at  small  emission 
angles.  But,  considering  the  whole  spectrum  as  made  up  from  a  'shear'  and  a 
'self'  noise  spectrum,  our  results  indicated  that  at  small  angles,  the  'shear' 
noise  spectrum  dominated  because  of  the  relatively  low  magnitude  of  the  'self' 
noise.  However,  as  the  angle  increased,  the  'shear'  noise  was  diminishing  as 
|cos0|.  At  right  angles  to  the  flow,  the  contribution  would  be  from  the  'self' 
noise  only  which  peaked  at  a  higher  frequency  than  the  ’shear’  noise  as  found 
experimentally  here. 

8.3  Acoustic  Power 


Finally,  an  estimate  of  the  total  acoustic  power  generated  by 
the  present  4  inch  jet  was  made.  Based  on  existing  knowledge  about  the  noise 
sources  distribution  in  a  jet,  an  equivalent  volume  of  turbulence  responsible 
for  the  noise  emission  was  obtained  whereby  the  total  noise  power  could  be 
obtained  by  multiplying  this  volume  by  the  noise  power  from  our  volume  ele¬ 
ment.  A  detailed  calculation  is  shown  in  Appendix  F.  The  result  is  plotted 
as  a  single  point  in  Fig.  40  along  with  a  replot  of  Fig.  7  of  Ref.  6.  Our 
noise  prediction  inferred  from  the  turbulence  measurements  is  seen  to  be  about 
an  order  of  magnitude  higher  than  the  extrapolation  of  the  curves  of  measured 
noise  obtained  by  others.  This  prediction  of  overall  jet  noise,  however,  is 
itself  a  tenuous  extrapolation  of  data  applying  to  a  single  volume  element. 

It  involves  a  number  of  assumptions  of  considerable  uncertainty,  and  it  should 
not  be  surprising  that  the  different  error  sources  might  accumulate  to  give  the 
observed  discrepancy. 
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IX.  CONCLUSION 


The  aerodynamic  noise  problem  first  solved  by  light hill  has 
been  reformulated  in  terms  of  its  spectral  characteristics  using  the  or.e- 
dimensional  Fourier  cosine  transform.  The  final  formalism  was  Just  as  re¬ 
vealing  as  other  workers'  results  using  tin  four- dimensional  Fourier  Transform. 
The  Doppler  effect  of  moving  eddies  resulted  automatically  from  the  mathematics 
and  the  convection  effect  could  be  distingui shed  as  a  combination  of  Doppler 
shift  and  amplification.  The  extent  to  which  differences  in  retarded  time 
could  be  neglected  was  explicit  in  this  formulation;  the  importance  of  retain¬ 
ing  the  differences  in  retarded  time  at  high  speed  was  demonstrated  following 
Ribner’s  work. 

The  main  objective  of  the  present  investigation  was  to  esti¬ 
mate  from  turbulence  measurements  the  characteristics  (i.e.,  intensity, 
directivity,  and  spectral  density)  of  the  noise  generated  by  a  unit  volume  of 
turbulence  in  the  mixing  region  of  a  model  jet.  To  accomplish  this,  two- 
point  space-time  correlations  of  both  the  turbulent  velocities  and  the  square 
of  these  velocities  have  been  obtained.  These  measurements  included  fine 
ordinary  u-component  and  also  the  components  at  4t°  and  60°  to  the  jet  axis. 

The  measured  quantities  approximate  the  actual  ir.teqrar.j 
the  Prcudman  form  of  Lighthill's  integral  connecting  turbulence  ard.  the  ••cicc 
it  radiates;  they  go  well  beyond  the  measurements  reported  els  own  ve  wh:  ..h 
rely  on  analogies  to  simulate  the  integrand .  The  requirement  to  deli  neat-*  a 
three-dimensional  integral  has  necessitated  a  great  volume  of  measuvouor4 s . 
Moreover  the  requirement  of  a  fourth  derivative  operation  on  the  experimental 
curves  has  put  a  premium  on  precision  of  measurement  and  data  processing. 

Where  comparable  measurements  exist  comparisons  have  bo  n 
made.  Both  the  intensity  of  the  u- component  and  the  eddy  convection  velocity 
agreed  well  with  those  measured  elsewhere,  Additionally,  comparable  space- 
time  correlations  are  qualitatively  similar  to  those  of  other  workers. 


Results  from  the  quantitative  estimate  of  the  noise  charactr- 
istics  suggested  that  the  basic  directivity  of  the  noise  generated  by  a  unL: 
volume  of  turbulence  varied  about  like  (l  +  1.6  jeose  j  ) .  The  absence  of  a 
valley  along  the  jet  axis  agreed  with  Ribncr's  finding  based  on  a  theoretic.' 
model;  it  disagreed  with  the  ' four- leaf  clover'  pattern  suggested  by  birht’.r! 
and  Lilley.  The  observed  dimple  in  experimental  jet  noise  patterns  along  tv 
jet  axis  is  due  to  another  effect,  namely,  refraction.  At  small  armies  zj 
the  jet  axis  the  'shear'  noise  inter.s.tty  was  found  to  be  more  than  *vi.e  the 
of  the  'self'  noise;  at  light  angles  the  'seif  noise  was  inferred  to  he  th.fr 
main  contributor.  The  'self  noise  spectrum  peaked  at  a  higher  frequency 
than  the  'shear'  noise  as  expected  from  theory.  Also  the  peak  frequenv  of 
the  present  unit  volume  spectrum  agreed  well  with  extrapolated  value  from 
acoustic  measurements  of  other  investigators.  However,  the  estimated.  •>* of  al 
acoustic  power  was  about  an  order  of  magnitude  higher  than  the  ext rap oj ated 
value  from  acoustic  measurements  based  on  the  U°  lav. 


!.  Total 
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In  addition,  the  present  investigation  had  demonstrate;:..  ‘-he 
usefulness  of  considering  the  'shear'  and  'self  noise  spearately  sorv-rhat 
in  the  manner  suggested  by  Ribner.  Looking  at  the  relative  importance  -f 
these  two  types  of  noise  at  different  angles  to  the  jet  axis,  one  can  tee  tic 
at  least  qualitatively  between  4 6°  to  $Kj°.  'The  acoustic  'paradox*  thau  the 


overall  noise  spectrum  peaks  at  higher  frequencies  for  larger  angles  from  the 
jet  axis.  This  is  contrary  to  the  deduction  from  the  Doppler  shift  of  a  single 
spectrum,  considering  the  eddies  as  moving  sources. 

Finally,  it  is  thought  that,  compared  with  earlier  work,  the 
present  investigation  has  taken  a  much  closer  step  to  a  quantitative  connection 
between  measurements  of  turbulence  and  measurements  of  jet  noise. 
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FIG.  1  CO-ORDINATE  SYSTEM 
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FIG.  2  NOISE  SPECTRUM  RADIATED  BY  UNIT  VOLUME  OF  JET 
TURBULENCE:  CONVECTIVE  PEAK  AMPLIFICATION  AND 
’DOPPLER’  SHIFT  ARE  GIVEN  BY  FACTORS  ©*AND  0"' 
RESPECTIVELY,  TO  GIVE  POWER  AMPLIFICATION 
(ADOPTED  FROM  REF.  3) 
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FIG.  4  CLOSE-UP  OF  HOT-WIRE  PROBES 


STAGGERED-HEAD  TAPE  RECORDER  PLAYBACK 
ASSEMBLY  FOR  INTRODUCING  TIME  DELAY 
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FIG.  6  FREQUENCY  RESPONSE  OF  ELECTRONIC  CORRELATOR 
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FIG.  7  AUTOCORRELATION  OF  SINE  WAVES  (CARRIED 
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FIG.  10  BLOCK  DIAGRAM  OF  THE  VELOCITY  SQUARED 
CORRELATION  MEASUREMENT  SYSTEM 


SPACE  CORRELATIONS  OF  BOTH  THE  VELOCITY  AND  THE  VELOCITY  SQUARED  IN  THREE 
OR  1 IIOCONAL  AXES;  '7,/DM.  L,/.i  -().  1J 
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FIG.  13  LONGITUDINAL  SPACE-TIME  CORRELATION  OF  u  f,/D*4,  f./D«0.  5,  U./a  *0.  13 
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FIG.  14  LONGITUDINAL  SPACE-TIME  CORRELATION  OF  u 
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LONGITUDINAL  SPACE-TIME  CORRELATION  OF  u_„  ;  t, /D*4,  ?»/ D*0.  5,  U  /•  -0.  13 


LOM'.n  IHJiNAL  SPACE  -  TIME  CUHIU.I  Al’IDN  OF  u T./D- 


‘juaioiwaoo  uoud|9jjoo 


LONGITUDINAL  SPACE-TIME  CORRELATION  OF  u45*;  1,1  D‘4.  ?«/ D*0.  5,  l)  /a  *0.13 


FIG.  18  LONGITUDINAL  SPACE-TIME  COH RELATION  OF 


FIG.  19  TYPICAL,  THREE  DIMENSIONAL  DISPLAY  OF  LONGITUDINAL  SPACE-TIME 
CORRELATION  OF  TURBULENT  VELOCITY 


FIG.  20  FIXED  AND  MOVING  FRAME  SPECTRA;  If/ 
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FIG.  21  TRANSVERSE  SPACE  CORRELATIONS  OF  u  .IN  2-AXIS; 
7,/D=4,  7i/D=0.  5.  U./a  =0.13  ° 


FIG.  22  TRANSVERSE  SPACE  CORRELATIONS  OF  uQ.IN  3*AXIS; 
7,/D=4,  ?i/D=0.  5,  U./a  =0.  13 


FIG.  23  TRANSVERSE  SPACE  CORRELATIONS  OF  U...IN  2- AXIS; 
?, /D=4,  f,/D=0.  5.  U./a  =0.  13  . 


TRANSVERSE  SPACE  CORRELATIONS  OF  u,*.IN  3-AXIS; 


FIG.  25  TRANSVERSE  SPACE  CORRELATIONS  OF  ufin.IN  2-AXIS; 
7,/D=4,  7j/D=0.  5,  U./a  *0.13 


FIG.  26  TRANSVERSE  SPACE  CORRELATIONS  OF  uRn.IN  3-AXIS; 
7)  /D=4,  7,/D=0.  5,  u./a  =0.13  b0 
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FIG.  28  TYPICAL  LONGITUDINAL  SPACE-TIME  CORRELATION  OF  TURNULENT  VELOCITY  SQUARED 
AFTER  MOVING-FRAME  TRANSFORMATION  USING  DATA  OF  FIG.  18 
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FIG.  29  FOR  THE  SHEAR  NOISE  CASES 

(a)  LEAST  SQUARE  CURVE  FITTING  TO  THE  INTEGRATED 
MOVING- FRAME  CORRELATION  FUNCTIONS 

(b)  THE  CORRESPONDING  FOURTH  DERIVATIVES 


FIG.  30  FOR  THE  SELF  NOISE  CASES 

(•)  LEAST  SQUARE  CURVE  FITTING  TO  THE  INTEGRATED 
MOVING- FRAME  CORRELATION  FUNCTIONS 
(b)  THE  CORRESPONDING  FOURTH  DERIVATIVES 


4th  Derivative  4th  Derivative  at  t  = 


FIG.  31  (a)  VARIATION  OF  FOURTH  DERIVATIVE  AT  t*0  OF  THE  FITTING 

POLYNOMIAL  (BY  LEAST  SQUARE)  VO  A  GAUSSIAN  FUNCTION 
WITH  THE  NUMBER  OF  DEGREE  OF  THE  POLYNOMIAL 
(b)  COMPARISON  BETWEEN  THEORETICAL  AND  COMPUTED 
FOURTH  DERIVATIVES 


Spectrum  (obtained  from  Fourier  cosine  transform ,  linear  scale) 


Frequency  -  Hz 


FIG.  32  COMPARISON  BETWEEN  THE  FOURIER  COSINE  TRANSFORM 
OF  THE  FOURTH  DERIVATIVES  OF  A  GAUSSIAN  FUNCTION 
AND  THE  APPROXIMATING  FUNCTION 
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FIG.  34  TYPICAL  COMPARISON  OF  THE  FOURTH  DERIVATIVES 
FROM  THE  THREE  TYPES  OF  CURVE  FITTING 
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FIG.  35  REPLOT  ON  A  LOG  LOG  SCALE  OF  THE  TAIL-  END  DATA 
OF  FIG.  29(a)  (NORMALIZED  BY  VALUE  AT  T  =0  ) 


FIG.  36  BASIC  DIRECTIVITY  OF  JET  NOISE  FROM 
TURBULENCE  MEASUREMENTS 
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(Q)  0*  CASE,  SHEAR  NOISE 


(b)  O' CASE,  SELF  NOISE 


FIG.  37  APPROXIMATING  THE  FOURTH  DERIVATIVE  CURVES 
BY  ANALYTICAL  FUNCTIONS 


Spectrum  (relative  intensity,  linear  scale) 


(b)  0*  CASE,  SELF  NOISE 

FIG.  38  ESTIMATED  SPECTRA  OF  SHEAR  AND  SELF  NOISE  RADIATED 
FROM  A'SLICE'OF  JET  AT  t/D*4,  0°  CASE 


LOG  ( fp#a|,)  SOUND  PRESSURE  in  db/cycl* 


LOG  (U) 

FIG.  39  REPLOT  OF  LEE'S  RESULTS  (  REF.  43  )  TAKEN  AT  25° 
FROM  JET  AXIS 
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ACOUSTIC  MEASUREMENTS  (  REPLOT  OF  FIG.  7,  REF.  6  ) 


APPENDIX  A 


Sensitivity  of  Inclined  Hot-Wire 

There  is  considerable  disagreement  as  to  the  directional 
sensitivity  of  a  hot-wire  (Reft  48,  49,  50,  51>  52,  and  53  to  quote  a  few). 

Of  these,  the  first  two  reported  the  usual  sine-law  response  for  angles 
greater  than  20°,  whereas  the  rest  reported  some  effects  of  the  tangential 
component.  A  complete  systematic  investigation  will  be  beyond  the  scope  of 
the  present  experiment.  Only  one  test  has  been  performed  with  parameters 
nearly  identical  to  those  that  would  be  used  in  subsequent  experiment,  i.e., 
linear  operation  mode  with  tungsten  wire  having  length- to-diameter  ratio  of 
230  and  overheating  ratio  of  0.3. 

The  test  was  performed  in  the  laminar  core  of  the  4  inch  jet. 
The  hot-wire  probe  was  mounted  in  a  rotatable  disk  assembly  as  shown  in  Fig. 
Al.  A  telescope  was  used  to  assure  that  the  hot-wire  was  at  the  center  of 
rotation.  The  inclination  of  the  wire  was  then  determined  from  the  angle  of 
rotation  of  the  disk.  The  initial  zero  position  was  determined  by  obtaining 
the  same  bridge  current  with  the  wire  inclined  at  a  positive  and  a  negative 
angle  to  an  arbitrary  fixed  datum.  The  mid  point  between  these  two  angles 
would  be  the  zero  datum. 

For  each  angle  of  inclination,  the  speed  (measured  with  a  pito 
static  probe)  of  the  jet  was  changed  and  the  bridge  current  versus  speed  was 
plotted.  Since  the  linear  operation  mode  was  used,  a  least  square  fit  was 
used  to  obtain  the  best  possible  straight  line.  Before  changing  to  a  differ¬ 
ent  inclination,  the  wire  was  first  brought  back  to  the  zero  position  and  the 
bridge  current  was  checked  for  the  same  maximum  speed.  From  these  bridge 
current  versus  speed  plots  for  different  inclination  angles,  the  sensitivity 
of  the  inclined  wire  for  three  representative  speed  ranges  was  obtained  and 
is  shown  in  Fig.  A2.  Results  indicated  that  there  was  no  systematic  devia¬ 
tion  from  the  sine  law  for  our  particular  case. 
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A2  SINE  LAW  RESPONSE  OF  INCLINED  CONSTANT- TEMPERATURE 
HOT-WIRE 


APPENDIX  B 


Effect  of  Low-Frequency  Cut-Off  on  Correlation  Functions 

At  tne  time  when  the  present  author  was  looking  into  this  j 
blem  he  was  unaware  of  similar  work  done  by  Wooldridge  and  Willmarth  (Ref. 
and  Bull  et  al  (Ref.  55)*  However  the  present  approach  takes  a  different 
form  and  is  based  on  a  relationship  between  a  filtered  and  a  non-filtered 
space-time  correlation  derived  by  Favre  et  al  (Ref.  56)  •  Moreover  the  eff 
on  the  shapes  of  the  space-time  correlations  are  considered  in  greater  det 

Suppose  the  turbulent  field  is  stationary  in  time,  the  spac 
time  correlation  between  the  fluctuating  velocities  at  two  points  separate 
by  a  distance  can  be  written  as 

R(1,t)  ■  UU7  (|,t)  =  J  |[u  +  U'(£,t)]*  -  [u  -  U'Q,t)]2| 


(B.l) 


The  filtered  space-time  correlation  is 


Rf(i>T)  ■  ufuf(i»T)  =  5  |[uf  +  -  [uf  -  u^(_|,t)-]2  j  (B.2) 

Now  if  the  frequency  response  (or  energy  response  to  be  more  exact)  of  the 
band-pass  filter  is  A(f),  then  by  definition 


>f  1  (|,t)]2  .  f  A(f)  0±(f)  df 

00 

[u  +  u*(1,t)]2  =  J  $>±(f)  df 


(B.3) 


(B.4) 


where  $t(f)  is  the  respective  power  spectral  density  function.  Let  us  de 
fine  the  following  autocorrelation  functions, 


R+  (t’>  =  [u  +  u'(£,t)]  [u(t')  i  U'  (_|,T  +  T')] 


(B.5) 


Then  using  the  Wiener-Khinchin  relationship 


4>+(f)  =4  f  R+  (t')  cos  27rfT’  dT’ 

J0  " 


Now  expanding  the  autocorrelation  functions,  we  have 

14  (t  1 )  =  R(0  ,t  '  )  +  R '  (0,t  '  )  i  [R(i,  T+  T*)  +  R(_|,t  - 


(B.6; 


(B.7! 


B1 


From  Eqs.  B.2  and  B.3, 


Vi’T)  =H  f  A(f)  *+(f)df  -  f  A(f)<J>  (f)dfl 
LJ0  J0  J 

*  f  A(f)df  f  R+(T')cos27rfT'dT' 

J°oo  J°oo 

-  [  A(f)df  J  R_(t')  cos27rfT'dT' 

J  0  Jo 

Inserting  R+(t’)  we  obtain 

00  00 

H  (i,r)  -  2  J  A(f)df  j  [R(1,t+t')  +  R(1,i-t’)]  cos27rfT'dT' 
The  above  equation  is  the  result  of  Favre  et  al  (Ref.  56). 


(B.8) 


(B.9) 


Since  we  are  interested  in  the  effect  of  low-frequency  cut¬ 
off,  we  assume  a  hypothetical  energy  response  of  the  system  to  be 


A(f)  =  (1  -  e 


2  /-  2 


■*7f, 


) 


(B.10) 


where  fc  is  the  cut-off  frequency  corresponding  to  a  2  db  fall  off.  This 
assumed  form  is  more  realistic  than  a  sharp  vertical  cut-off.  Inserting 
(B.10)  into  (B.9), 

00 

Rf(!>T)  =  2  JJ  [R(I»  T+  T’)  +  r(1>t  -  t’)]  cos27^fT,dfdT, 

0 

-  2 J°°  p(i,T+T')  +  R(],t.t')]J  e"f2/fc  cos27rfT,dfdT'  (B.ll) 

0 

The  first  term  is  just  the  unfiltered  space-time  correlation,  as  can  be  seen 
from  (B.9)  with  A(f)  =  1.  Therefore, 

Rf(i,T)  =  R(I,t)  -  2  J  e'f  /fc  dfj  [R(5,t+t')  +  R(1,t-t')]  cos27TfT'dT' 

(B.12) 

Note  that  even  when  t  =  0  the  space  correlation  will  also  be  affected  by  the 
low-frequency  cut-off.  The  above  equation  is  useful  in  estimating  the  filter¬ 
ed  correlation  from  the  non- filtered  one,  but  the  reverse  procedure  would 
involve  solving  an  integral  equation.  By  a  change  of  variables,  one  can  show 
that 


J'  [r(^,t+t ' )  +R(_I,t-t')]  cos277"fT  ' dT ' 

0 

00 

=  J  R(5_,T)  [cos27rfTcos27rfT  +  sin27nfTsin277fT] 


dT 


(B.13) 


—  00 


BP 


Therefore 

=  R(i>T) 


2 


[cos277fTcos27rfT  +  sin27rfT 


x  sin27rfr  ]dT 


(B.14) 


Now  one  can  identify  the  two  terms  in  the  last  integral  as  the 
real  and  imaginary  parts  of  the  cross-spectral  density  associated  with  R(_|,t), 
i.e.  „ 


Hf(i,T)  =  R(1,t)  -  2  J 

0 

=  R(|,t)  -  2  J” 


where 


_f2 

e  C  df(27T  ^(J^uJcosut  -  277<>j(_^,oj)  sinurr) 
-oj2/u§ 


dw  [  ,td)  |  COS  a  COSWT 


=  R(£,t)  -  2  e 


f 


l*(I  ,w)  |  sin  a  sinurr  ] 
-u)2yw2 


|<t>(_|,u)  |  cos(uvr  +  a)  dw 


|2  =  4>r2(£,u)  +  <t|(£,u) 


a  =  tan 


R 

4>I 


-1 


4V, 


(B.15) 


If  we  follow  Bull  et  al's  assumption  that  the  cut-off  is 
sharp  and  |<t>(_£,uj> }  is  constant  from  the  cut-off  frequency  down  to  zero  fre¬ 
quency,  we  obtain  instead 


>>c 

Rf(6,T)  =  R(£,t)  -  2  i4>(l,w)|0  /  cos(wt  +  a)du  (B.l6) 

Defining  a  convection  velocity  Uc  by  a  =  -t]w/Uc  and  assuming  it  to  be  in¬ 
dependent  of  frequency,  the  result  becomes 


Rf(i,T)  =  R(£,T)  -  2  |4>(£,u)l0 


sin  uc  (t  -  tt~  ) 
uc 


(t  - 


(B.17) 


which  is  identical  to  that  obtained  by  Bull  et  al.  Although  the  convection 
velocity  defined  here  is  slightly  different  from  that  defined  in  the  text, 
they  can  be  considered  as  equivalent  to  a  first  approximation.  If  this  is 
the  case,  then  on  transforming  to  the  moving  frame  the  error  term  in  (B.17) 
will  be  independent  of  t  as  Uct  -  ^  becomes  -§mi.  It  follows  that  the  lov- 
frequency  cut-off  will  not  affect  the  fourth  time-derivative  of  the  moving- 
frame  space-time  correlation.  However,  it  will  affect  our  length  scales 
especially  the  longitudinal  one. 


Rf(i,0)  =  R(i,0)  -  2 


From  Eq.  B.l4,  when  t  =  0  we  obtain 

.  r  .-**/**  --  r" 


'  f  R(i.i) 

kLoo 


cos27rfT  dT 


(B.l8) 


Let  us  first  consider  the  transverse  correlations.  A  plausible  form  for  the 
space-time  correlation  will  be 

.m2 


R(iitT)  *  R(81,0)  e~ai 
Then  it  can  be  shown  fairly  easily  that 

RfU^O)  =  R(ii,0)  (l  - 


i  i  -  2,3 


(B.19) 


(B.20) 


The  above  equation  indicates  that  the  low-frequency  cut-off  will  affect  the 
magnitudes  of  the  transverse  space  correlations  but  not  their  characteristics. 
In  other  words,  if  the  transverse  space  correlation  is  non-negative,  it  will 
remain  non-negative.  The  situation  is  quite  different  for  the  longitudinal 
space  correlation  as  seen  in  the  next  paragraph. 

In  order  to  simplify  the  mathematics,  we  have  chosen  a  'frozen' 
convective  form  of  the  space-time  correlation  for  the  longitudinal  case,  i.e., 


r(51,T)  =  e 


-A(5i  -  UCT)2 


(B.21) 


H(81,0)  =  e  1 

With  this  assumed  form,  Eq.  B.l8  can  be  integrated  to  give 


*1  1  ' 


RfU-^O)  =  R(?1,0)  1  - 


_  a 
2  + 


(B.22) 


One  can  see  that  for  large  enough  the  bracketed  term  can  become  negative. 

In  other  words,  even  if  R(l1,0)  is  non-negative  as  in  the  assumed  form  here, 
the  low-frequency  cut-off  will  render  it  negative  at  large  1^*  In  fact  the 
resultant  space  correlation  Rf(li,0)  will  have  equal  positive  and  negative  areas 
This  can  be  shown  by  integrating  Eq.  B.22  to  infinity  with  respect  to  . 

Thus  the  ordinary  definition  of  length  scale  defined  by  this  integral  will 
no  longer  be  meaningful  because  it  will  be  zero.  Thus  one  must  be  careful 
in  interpreting  the  length  scales  obtained  by  different  workers:  the  trunca¬ 
tion  of  the  integral  in  conjunction  with  low-frequency  cut-off  effect  will 
lead  to  a  variety  of  "scales". 

Since  we  have  used  the  'frozen'  convected  pattern,  the  remarks 
of  the  last  paragraph  can  also  be  applied  to  autocorrelations  and  time  scales. 
This  is  possible  because  there  is  interchangibility  between  autocorrelation 
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and  longitudinal  space  correlation  in  a  'frozen'  convected  turbulence  field. 

Except  for  the  case  of  autocorrelation,  it  is  fairly  difficult 
to  correct  for  the  low-frequency  cut-off  effect  because  it  will  involve  mea¬ 
surements  of  the  cross-spectral  density  function  as  indicated  by  Eq.  B.15.  In 
the  forthcoming  Reference  57,  an  actual  case  is  presented  in  which  an  auto¬ 
correlation  function  has  been  corrected  for  the  low-frequency  cut-off  from  a 
knowledge  of  the  spectral  density  function.  Here  the  best  we  can  do  is  to 
obtain  an  estimate  of  its  effect  on  the  length  scales. 

By  approximating  our  data  by  the  assumed  forms  of  the  space- 
time  correlations,  it  was  estimated  that  our  two  transverse  length  scales 
might  be  off  by  about  10 $  and  that  the  longitudinal  length  scale  might  be 
lower  by  as  much  as  20$  if  we  limited  the  scale  integral  to  the  positive  area 
only. 
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APPENDIX  C 


Functions  Used  in  Least  Square  Curve  Fitting 
The  functions  used  are  of  the  form 

A  e“(ao  +  a2  T2  +  &4  T4  +  ...) 


where 


T  =  300. Oxt  second 


'shear'  noise  cases 


0°  Case 

45°  Case 

60°  Case 

A 

0.766 

0.674 

0.576 

ao 

0.130350E-01* 

0.157922E-01 

0.238837E-01 

a2 

0. 11072 3E +01 

0.123883E+O1 

0.148080E+01 

a4 

-0.279533E+01 

-0.327688E+01 

-0.412232E+01 

0.511338E+O1 

0.566124E+01 

0.800265E+01 

a8 

-0.483750E+01 

-0.510798E+01 

-0.820410E+01 

a10 

0.222496E+01 

0.226621E+01 

0.413337E+01 

a12 

-0.39^520E+00 

-0.390888E+00 

-0.806750E+00 

'self' 

noise  cases 

0°  Case 

45°  Case 

60°  Case 

A 

0.408 

0.320 

0.279 

ao 

0.159178E-01 

0.188755E-01 

0.304724E-01 

a2 

0.848186E+01 

0.664074E+01 

0.185413E+02 

a4 

-0.155508E+03 

-0.800311E+02 

-0.353693E+03 

a6 

0.164690E+04 

0.859596E+O3 

0.277588E+04 

a8 

-0.762908E+04 

-0.497161E+O4 

0.111248E+03 

a10 

0.770442E+04 

0.957882E+04 

-0.836756E+05 

a12 

0.516812E+05 

O.276362E+05 

0.125948E+06 

al6 

-0.167787E+06 

-0.142374E+06 

0.127532E+07 

0.146786E+06 

0.161550E+06 

-0.343257E+07 

-1 


*  0.130350E-01  =  0.130350x10 


APPENDIX  D 


Differences  in  Emission  from  Unit  Volumes  Around  the  Jet  at  Same  Radial  Distance 

from  Jet  Axis 


Although  there  is  axis- symmetry  in  the  overall  jet  noise  pattern, 
the  noise  radiated  in  a  certain  direction  from  volume  elements  around  the  jet 
at  a  given  downstream  position  will  not  be  the  same.  In  an  unpublished  work, 
Ribner  has  taken  into  account  the  dependence  of  radiation  from  volume  elements 
on  this  azimuth  angle  0  (see  figure  below) 


Looking  down  on  the  jet  exit. 
Observer  on  1-2  plane  making  9  with 
1-axis. 


and  obtained  the  following  corrected  theoretical  directivity 


a3i(x,0) 

" 


1  +  cos^0  +  cos20  sin20  sin2! 


As  0  varies  from  0  to  2i r  ,  this  varies  sinusoidally  between 

Jl  O 

1  +  cos  0  and  1  +  cos^O 
so  that  the  0-average  is 


d3l(X,0) 

<  - > 

av.0 


i  + 


cos^O 


^0  +  cos^9 


Thus  it  is  seen  that  for  the  theoretical  model  the  overall 
directivity  does  not  differ  markedly  from  his  original  (l  +  cos*4,  )  pattern 
taken  at  0  =  0  (Refs.  3  and  13) •  In  the  absence  of  other  information  we  shall 
assume  a  similar  relationship  for  the  experimental  resuits  herein;  that  is,' 
(with  other  assumptions),  the  directional  pattern  emitted  from  a  volume  ele¬ 
ment  at  ,  as  inferred  for  the  plane  containing  _]}  and  the  jet  axis,  is 
essentially  the  directional  pattern  inferred  for  the  jet  as  a  whole. 
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APPENDIX  E 


Spectrum  Functions  Obtained  from  Truncated  Autocorrelation  Functions 


One  of  the  difficulties  in  experimental  work  with  random  fluctua¬ 
tions  is  that  it  is  impossible  to  obtain  a  complete  form  of  the  autocorrelatior 
function  up  to  infinite  time  delay;  the  function  is  always  'truncated'  at  a 
finite  time  delay.  If  one  is  interested  in  the  spectrum  function  obtained 
from  the  autocorrelation  by  the  Fourier  cosine  transform,  some  idea  of  how 
this  truncation  of  the  autocorrelation  function  affects  the  final  spectrum 
function  should  be  known. 


A  hypothetical  case  is  a  pure  cosine  function.  The  corresponding 
spectrum  function  will  be  a  delta  function  centered  at  the  same  frequency.  If 
this  cosine  function  terminates  at  a  finite  time  t  ,  the  corresponding  spec¬ 
trum  function  will  be  of  the  form  sin(w-oo0)T/2(u)-<J0)  which  resembles  a  damped 
oscillation  centered  at  tjQ.  This  function  has  negative  values  and  approaches 
the  delta  function  for  large  t. 

A  more  realistic  autocoi relation  function  is  the  exponential 
cosine  function  which  has  both  the  infinite  and  finite  Fourier  cosine  trans¬ 
forms.  Figure  B1  again  shows  the  occurance  of  negative  values  and  a 
general  broading  of  the  spectrum  function  for  the  truncated  autocorrelation 
function.  However,  the  position  of  the  peak  frequency  is  not  affected. 
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FIG.  El  EFFECT  OF  TRUNCATED  AUTOCORRELATION  FUNCTION 
ON  SPECTRUM 


APPENDIX  F 


Estimation  of  Acoustic  Power 


Based  on  a  theoretical  model,  Ribner  (Ref. 3)  has  shown  that  the 
acoustic  intensity  emitted  in  direction  0  from  a  volume  element  is  (as  corrected 
in  Appendix  D) 


d3l(X) 


*3n 


|X|! 


a 

o 


( 


2 

cos  9 


cos  0 


■) 


(F.l) 


which  can  be  re-written  as 


Values  with  a  designate  values  at  arbitrary  radial  position  T)^ 
whereas  values  without  a  refer  to  that  value  of  ^  where  33t(x)/^3jj  peaks^ 
for  given  The  frequency  Uf  is  a  typical  radian  frequency  in  the  turbu¬ 

lence;  it  is  defined  by  approximately  matching  a  function  exp(-wf/r/)  to  the 
experimental  autocorrelation  of  u1  in  the  convected  reference  frame. 


For  a  given  jet  at  a  given  operating  condition,  f  and  g  are 
functions  of  position  only:  f  =  r( ) ,  g  =  g(T)2jT)-,).  The  function  g(rj2,T]1) 
is  a  measure  of  the  radial  distribution  of  the  effective  sound  sources.  The 
integral  of  g(l2»r)l)(^  measures  the  effective  area  S(t]^)  occupied  by  the  sound 
emitters  in  a  'slice'  of  jet  dT)^  at  t;-^,  considered  as  replaced  by  emitters  of 
the  (radially)  uniform  strength  f(i]^).  Thus  we  can  write 


dl(x) 

dT)i 


2_+  cos^9+  cos^8 


3(1!)  (F.3) 

(peak  value  atT}!) 


Four  our  case,  we  have  measured  experimentally  the  value  on  the 
right-hand  side  in  the  more  exact  form  at  T]!  =  4d.  It  is  only  necessary  to 
estimate  3(1-^  from  the  g  function  at  the  same  location.  The  effective  area 
will  be  an  annular  one  as  shown: 


FI 


Now  across  the  mixing  region  only  uj_  and  will  vary 
ippreciably.  The  variation  of  was  obtained  from  the  present  experiment 
whereas  the  variation  of  had  to  be  inferred  from  Ref.  8,  assuming  simil- 

irity  between  the  two  Jets.  The  effective  area  was  estimated  to  be 

8(1^)  -  0.12  tj-l  D  (F.4) 

(  *  0.48  for  our  case) 

/hich  agreed  with  the  result  of  Ref.  12  using  a  slightly  different  function. 

The  power  generated  by  a  slice  of  Jet  was  then  obtained  by  inte- 
jrating  the  intensities  at  different  directions  over  the  surface  of  a  sphere. 
!he  experimental  basic  directivity  of  (l  +  2.6/cos0/)  was  used  and  the  con¬ 
nection  factor  C”5  was  replaced  by  0~5  for  convenience  in  integration. 

Since  'slices'  of  jet  in  the  mixing  region  emit  the  same  noise 
rawer  and  the  emission  falls  off  very  abruptly  (about  like  t^-?)  beyond  six 
liameters  or  sc  (Refs.  12,  14,  and  15),  the  total  acoustic  power  was  roughly 
>btained  by  multiplying  the  result  for  our  unit  length  'slice'  by  a  length 
if  six  diameters.  The  estimated  power  generated  by  our  jet  was  about  1.46  x 
.0"5  watts  according  to  this  estimation. 
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